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Assembly  and  Application  of  a  Shock  Data  Analysis  System 


1.  Introduction 

Measured  data  from  shipboard  shock  environments  usually  consist  of 
analogue  magnetic-tape  recordings  of  acceleration  and  strain  versus  time. 
Usually,  little  analysis  is  required  of  strain  recordings  save  spectral  analy¬ 
sis  and  perhaps  sums  and  differences.  From  the  accelerometer  recordings, 
however,  it  is  usually  necessary  to  extract  velocities  and  displacements,  and 
often  shock  spectra,  as  well.  Shock  data  analysis  facilities  tend  to  run  to 
two  extremes  of  size  and  capacity.  At  one  end  is  the  large  digital  processing 
facility  where  hundreds  of  records  (from  an  at-sea  ship  shock  test,  for 
example)  are  processed  in  batch,  usually  after  the  test  is  completed.  At  the 
under  end  is  the  small,  often  analogue,  facility  which  can  do  limited  analysis 
of  a  small  number  of  records,  although  perhaps  very  slowly.  The  latter  may  be 
used  to  advantage  for  analyzing  a  few  selected  critical  records  in  the  field, 
or  in  the  shock  simulation  lab. 

A  need  exists  for  a  capability  in  the  middle  ground  -  a  system  capable  of 
doing  analysis  on  a  modest  number  of  records  in  a  reasonable  time,  which 
accepts  analogue  magnetic  tapes  as  input  and  produces  report-ready  graphs  and 
tables  as  output,  and  which  is  field-transportable  if  necessary. 

This  report  describes  such  a  system  which  has  been  assembled  with 
standard,  off-the-shelf  instruments  which  are  marketed  nationwide  and  easily 
available.  It  is  not  in  any  way  optimized  for  the  specific  task  of  shock  data 
analysis,  hence  may  not  be  competitive  in  speed,  cost,  etc.  to  a  system  which 
has  been  optimized  in  hardware  and/or  software  to  this  specific  application. 

It  furnishes  an  example  of  what  can  be  done  along  these  lines  with  virtually 
no  expenditure  of  R  &  D  effort,  and  is  an  example  in  that  any  of  its  component 
instruments  could  be  replaced  by  a  different  supplier's  product  with  little 
change  in  operation  or  performance. 

2.  Background 

2.1.  Shock  motions  may  be  described  in  terms  of  acceleration 
velocity  or  displacement- time  histories,  and  in  all  three  parameters  may 
vary  greatly  in  waveform  complexity.  Simpler  shock  motions  may  be  closely 
defined  by  a  few  waveform  parameters,  and  these  may  also  provide  a  reliable 
basis  for  comparison  of  shock  motions.  Shipboard  shock,  however,  is  extremely 
complex.  Motions  are  typified  by  waveforms  of  high  frequency  range,  high 
dynamic  range  and  long  duration  which  cannot  be  defined  adequately  by  a  few 
parameters.  The  additional  descriptions  found  most  useful  for  a  motion  of 
this  type  are  the  shock  spectrum,  which  may  be  interpreted  as  a  measure  of  the 
effects  of  a  shock  motion,  and  the  design  shock  spectrum,  which  provides  the 
data  needed  to  design  a  structure  to  withstand  it.  Each  provides  a  basis  for 
comparing  shock  severities,  and  applies  equally  well  to  shock  motions  of  all 
degrees  of  complexity  of  waveform.  The  general  nature  of  shock  motions,  shock 
spectra  and  design  shock  spectra  are  discussed  in  Refs.  (1)  and  (2). 

2.2.  While  a  shock  motion  may  be  identified  equally  well  by  its  acceler- 
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ation,  velocity  or  displacement  waveform,  one  of  them  will  usually  be  most 
convenient  for  a  given  application.  However,  in  most  situations  only  acceler¬ 
ation  can  be  measured  properly.  This  is  almost  always  the  situation  in  the 
case  of  shipboard  shock,  and  the  bulk  of  data  in  this  field  will  consist  of 
records  of  acceleration  and  strain  versus  time,  with  perhaps  a  small  fraction 
being  records  of  velocity  versus  time.  Thus,  an  analysis  system  for  shipboard 
shock  motion  records  must  be  able  to  produce  time  histories  of  acceleration, 
velocity  and  displacement,  as  well  as  shock  spectra  and  design  shock  spectra, 
given  input  signals  representing  acceleration  -  or  velocity-time  histories, 
and  preferably  from  displacement-time  histories  as  well.  Since  shipboard 
shock  motions  are  among  the  most  demanding  in  terms  of  functional 
requirements,  an  analysis  system  capable  of  handling  them  will  also  be  able  to 
handle  shock  motions  of  most  other  types. 

3.  Analysis  System  Characteristics 

3.1.  The  principal  application  of  the  system  is  the  analysis  of  shock 
motions  measured  on  laboratory  shock  machines.  This  includes  those  which 
generate  simple-pulse  shock,  but  most  prominent  are  those  which  generate 
shipboard  shock  type  environments:  the  Lightweight  Shock  Machine  (LWSM),  the 
Mediumweight  Shock  Machine  (MWSM) ,  the  Floating  Shock  Platform  (FSP)  and  the 
Large  Floating  Shock  Platform  (LFSP)  (Ref.  3).  Tests  on  these  machines  will 
be  lightly  instrumented  compared  to  at-sea  tests  of  ships  and  submarines,  or 
some  of  the  more  specialized  test  vehicles.  The  typical  instrumentation  suite 
may  be  a  few  channels  on  the  simpler-pulse  machines  and  the  LWSM,  up  to  a 
dozen  on  the  MWSM  and  FSP  and  a  few  dozen  on  the  LFSP,  while  on  some  occasions 
several  times  as  many  channels  may  be  required.  Normally,  about  one- third  of 
the  total  channels  will  fall  in  the  category  of  equipment  input  motions,  while 
most  of  the  remainder  will  be  equipment  response  motions  and  structural 
strains. 

3.2.  The  function  of  the  analysis  system  is  to  accept  analogue  input 
signals,  perform  whatever  processing  is  required,  and  produce  output  represen¬ 
tations  of  the  desired  quantities,  scaled  in  the  appropriate  engineering 
units,  in  graphical  or  tabular  form,  and  in  video  or  hard-copy  display.  All 
signals  require  A-D  conversion  and  storage,  scaling  and  display.  All  motion 
signals  require  in  addition  computation  and  display  of  acceleration,  velocity 
and  displacement,  each  properly  scaled,  which  requires  the  functions  of 
integration  and  possibly  differentiation.  Input  motion  signals,  and  some 
response  motion  signals,  also  require  the  computation  of  raaxiraax  and  residual 
shock  spectra.  This  last  computation  is  by  far  the  most  complex  and  time- 
consuming,  but  is  needed  for  only  a  few  channels.  For  a  suite  of  a  dozen 
channels  of  assorted  information,  the  system  will  perform  a  complete  analysis 
in  an  hour  or  two. 


Shock  spectrum  analyzers  are  available  commercially,  and  are  typically 
compact  units  and  very  fast.  As  stand-alone  systems,  they  are  lacking  in  that 
they  cannot  do  a  complete  analysis  (integration  and  scaling,  etc.)  and  have  no 
display.  As  components  for  calculating  shock  spectra  they  are  wanting  in  that 
they  cannot  compute  undamped  spectra  and  their  analysis  frequencies  are  fixed. 


3.3.  Since  the  analysis  system  is  an  assembly  of  component  units  designed 
for  use  as  independent,  general  purpose  instruments,  it  has  a  great  degree  of 
flexibility.  With  different  switch  settings,  it  would  function  equally  well 
for  analyzing  earthquake  motions  or  blast  loadings.  A  brief  description  of 
each  functional  block  of  the  system  and  its  capabilities  follows. 

3.3.1.  Magnetic-Tape  Recorder 

The  primitive  form  of  shock  data  records  is  usually  FM 
magnetic  tape.  The  recorder  used  with  the  analysis  system  has  14- 
track  record/reproduce  capability,  plus  edge  tracks  for  voice  annota¬ 
tion.  The  normal  configuration  has  two  housekeeping  channels,  one  of 
them  DR  for  a  tape-speed  servo  reference  signal  and  the  other  an  FM 
channel  devoted  to  time-code  in  IRIG  B  format,  and  12  FM  data 
channels.  All  FM  channels  are  IRIG  Intermediate  Band  I.  Normal 
recording  speed  is  30  ips,  giving  a  nominal  bandwidth  of  0  -  10 
KHz.  The  recorder  can  be  operated  at  tapes  speeds  from  15/16  to  120 
ips,  allowing  the  traditional  analysis  tricks  of  data  speed-up  and 
slow-down;  record  density  can  be  changed  by  rearranging  plug-in 
components  on  the  electronics  cards.  Output  level  is  normally  set 
for  lv_  rms  for  ±40%  peak-to-peak  (sinusoidal)  modulation.  This 
particular  unit  was  selected  because  it  performs  a  set  of  canned 
calibration  and  set-up  routines  under  control  of  a  built-in  micropro¬ 
cessor,  and  because  it  is  a  portable  unit.  Both  features  are 
convenient  for  field  operations. 

3.3.2.  Auxiliary  Equipment 

The  principal  auxiliary  equipment  units  are  a  time-code  gener¬ 
ator/reader,  a  direct-writing  oscillograph,  and  a  storage  CRD.  In 
generator  mode,  the  time-code  unit  produces  an  IRIG  B  code  for 
recording  in  the  magnetic  tape.  In  reader  mode,  it  monitors  the 
recorded  time-code  and  produces  a  pulse  when  a  chosen  time  value 
occurs.  This  time  can  be  selected  to  the  millisecond  and  provides  a 
start  signal  to  the  transient-storage  unit  (discussed  below)  and 
facilitates  synchronization  of  the  recorded  data  signals.  The  other 
units,  the  CRO  and  oscillograph,  provide  quick-look  disDlay  of  the 
recorded  data  signals  and  time-code.  The  CRO  is  a  two-channel, 
single  time-base  type  similar  to  many  currently  available.  The 
oscillograph  produces  strip  charts  of  up  to  12  channels  with  a  band¬ 
width  of  0  -  5  KHz  at  paper  speeds  from  0.1  to  160  ips. 

3.3.3.  Data  Filter 

The  data  signals  output  from  the  tape  recorder  are  passed 
through  an  8-channel  data  filter.  Each  channel  has  a  low-pass 
characteristic  with  cutoff  frequency  selectable  from  1  Hz  to  15  KHz, 
and  is  direct-coupled.  The  unique  features  of  this  unit  are  that  the 
filter  characteristic  is  optimized  for  linear  phase-shift  within 
pass-band,  minimizing  envelope  distortion,  and  that  phase-shift 
between  channels  is  closely  controlled  at  no  more  than  2°  between  any 
pair.  The  filter  characteristic  also  features  a  cutoff  rate  of  43 
db/octave,  so  that  aliasing  is  negligible  for  sampling  frequencies  of 
4,2  times  the  cutoff  frequency  or  higher.  Each  channel  has  unity 
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gain,  and  can  accept  inputs  of  up  to  10  _V  «  As  used  in  the  analysis 
system,  the  cutoff  frequency  is  normally  set  to  lKHz. 

3.3*4.  Transient  Storage  Unit 

The  transient-storage  unit  contains  seven  independent  AD  and 
digital  storage  channels  with  a  common  tine-base.  Each  channel  can 
be  operated  in  an  inhibit  mode  which  preserves  the  existing  stored 
signal.  The  unit  can  be  operated  in  a  self-triggered  mode,  where  the 
signal  in  any  channel  can  be  used  to  initiate  the  capture  of  the 
signals  in  all  channels  which  are  not  inhibited.  As  used  in  the 
analysis  system,  the  capture-initiate  signal  is  derived  from  the 
time-code  reader  to  preserve  commonality  of  time-base  between  succes¬ 
sive  batches  of  signals.  Upon  initiation,  each  channel  captures  4096 
10-bit  samples  of  its  signal,  for  a  dynamic  range  of  60  db.  As  used 
in  the  analysis  system,  the  channel  full-scale  range  is  ±  2  V  ;  it 
can  be  varied  from  ±0. 1  to  ±5V^.  The  sampling  period  used  is  100 
microsec.,  for  a  total  record  length  of  400  millisec.;  the  sampling 
period  can  be  varied  from  5  microsec.  to  2.5  sec.  Each  channel 
provides  outputs  of  the  digitized  "staircase"  and  smoothed  analogue 
versions  of  the  stored  signal  waveform.  As  used  in  the  analysis 
system,  the  signal  stored  in  any  selected  channel  is  read  out  on  the 
IEEE-488  bus  in  bit-parallel,  character-serial  ASCII  code.  This 
particular  unit  was  selected  because  of  its  7-channel  capacity  and 
the  fact  that  each  A-D  converter  and  storage  channel  is  functionally 
independent  from  the  others. 

3.3.5.  Processor  and  Display  Unit 

The  processor  and  display  unit  is  a  large-screen  desk-top 
computer  and  peripherals.  The  computer  has  64K-byte  memory  and  a 
resident  BASIC  interpreter,  leaving  some  53K-bytes  of  usable 
memory.  Additional  ROM  packs  provide  signal  processing  routines, 
including  Integration  and  differentiation,  and  text-editing 
capability.  The  principal  display  is  a  built-in  large  screen  CRT. 

The  peripherals  include  a  hard-copy  unit  which  gives  a  permanent 
record  of  the  CRT  display  with  reduced  size,  and  a  digital  X-Y 
plotter  which  gives  high  quality  graphs  up  to  11  x  17  in.  Other 
peripherals  are  a  high  capacity  disc-file  manager,  which  allows  a 
large  library  of  programs  to  be  resident  on  floppy  disc,  and  an 
additional  tape-cassette  drive-unit,  allowing  more  flexibility  in 
data  transfer.  All  of  these  components  are  produced  by  the  same 
manufacturer,  and  are  configured  and  interfaced  to  form  a  compatible 
system.  This  (processor  and  display)  unit  was  selected  because  it 
possesses  the  requisite  computational  power  and  has  an  excellent 
repertoire  of  graphics  routines.  Other  features  are  good  man-machine 
interface,  BASIC  interpreter  and  large-screen  display.  These  seem  to 
assist  the  Intended  user  profile,  engineering/technical  personnel 
without  extensive  familarity  with  computing  systems,  in  learning  to 
operate  this  analysis  system  with  a  miniman  of  instruction.  The 
existence  of  a  comprehensive  assortment  of  proven  peripherals  was 
also  a  factor  of  consideration.  Another  feature  of  the  computer  is 
that  it  can  also  be  operated  as  a  high-speed  (9600  baud)  intelligent 
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terminal,  should  operation  with  a  larger  host  computer  be 
desirable.  The  last  peripheral  in  the  system  is  a  high-speed 
printing  terminal  from  another  manufacturer.  When  used  with  the 
ana  .system,  it  functions  as  a  moderately  fast  line-printer. 

3.4.  The  functional  block  diagram  of  the  analysis  system  is  in  shown  in 
Fig.  1. 

4.  Computation  of  Shock  Spectra 

4.1.  Recalling  that  the  shock  spectrum  of  a  shock  motion  may  be 
interpreted  as  the  maximum  relative  displacements  of  a  set  of  massless  single- 
degree-of-freedora  systems  whose  bases  are  forced  to  follow  the  motion,  it  is 
apparent  that  shock  spectra  may  be  derived  for  any  degree  of  damping.  For  the 
purpose  of  comparing  shock  motion  environments  and  specifying  design  shock 
spectra,  the  undamped  version  is  preferred,  since  the  damped  spectra  are  less 
sensitive  to  variations  in  the  input  motion.  Shock  spectra  can  be  computed 
from  the  differential  equation  of  motion,  or  from  the  corresponding  integral 
equation.  Either  form  can  be  expressed  to  involve  the  input  motion  in  terms 
of  acceleration,  displacement  or  velocity.  The  procedure  used  here  is  that  of 
numerical  solution  of  the  differential  equation  relating  relative  displacement 
to  input  accleration.  The  computational  procedure  requires  the  input  to  the 
program  as  the  acceleration-time  history  of  the  shock  motion.  The 
computational  procedure  uses  the  algorithm  described  in  Ref.  (4). 

The  equation  of  motion  for  an  undamped  oscillator  of  frequency  o  is 
( 1)  X  +  Ui*X  =  -Z 


with  the  solution  for  the  relative  displacement  X 
(2)  X(t)  =  XQ  cos  ait  +  sin  tot  -  /£z(t)  sin  u>(t  -  t)  dt . 

The  numerical  integration  equation  derived  by  G.  O'Hara  (Ref.  4)  is 


(3) 

Xn+^w  =  Xnu>  cos  u>h  4-  Xn  sin  o>h 


3n(  1  -  cos  u)h)  _  / 1  +  cos  uh 

u)h  5n-l'  2u)h 


sin  u)h\ 

STZTT^) 


with  a  parabolic  approximation  for  the  digitized  input  function 
here  Sn  "  Zn+1  -  Zn 


Z  (t); 


and  S*n-1 
(note  that  SZn_j 

so  that 


Sn  "  sn-l  *  zn+l  "  ^n  +  Vl 
is  the  second  difference,  not  a  square) 
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In  these  equations,  h  is  the  time  interval  of  integration. 

By  eliminating  Xn  in  (3),  R.  Bort's  implementation  of  equation  3  (Ref.  5) 
is 

(5)  X(t  +  h)  =  2  X(t)  cos  wh  -  X(t-h)  +  h  [Z  (t  +  h/2)  -  Z(t  -  h/2)] 

Equation  (5)  has  been  applied  in  this  program  because  it  consists  of  only 
three  additions.^  As  implemented  here,  the  term  in  brackets 
[Z  (t  +  h/2)  -  Z  (t  -  h/2)]  is  identified  as  hZ  (t) ,  a  sample  point  from  the 
digitized  input  acceleration  record.  Before  processing  input  data  the 
2  cos  u>h  term  is  determined  and  stored  in  array  C.  All  scaling  factors  are 
deferred  to  the  end  of  the  calculation  for  each  analysis  frequency. 

In  order  to  cut  down  calculation  time,  the  fact  that  lower  analysis 
frequencies  require  fewer  periods  to  reach  peak  response  than  higher 
frequencies  is  employed.  The  number  of  periods  per  analysis  frequency 
considered  in  the  processing  routine  is  given  by 

0  =  3"^og  T 

where  F  =*  oj/2it  . 

The  second  procedure  to  speed  up  the  program  selects  data  points  with  respect 
to  cutoff  frequency  and  analysis  frequency.  The  calculation  time  has  been 
reduced  this  way  from  about  50  min.  to  8  min.  while  producing  only  a  minor 
error  (see  Fig.  2A,B). 

5.  System  Operator's  Manual 

5.1.1.  It  is  assumed  that  a  time  record  of  the  shock  motion  - 
including  a  calibration  step  -  is  available  on  an  analog  tape  and  a 
time  code  is  recorded  on  a  separate  track  of  the  tape.  Playing  back 
the  record  (or  calibration  step)  and  the  time  code  on  the  Visicorder, 
determine  the  millisecond  preceding  the  starting  time  of  the 
record.  Set  the  time  code  generator  to  open  the  time  gate  and 
trigger  the  transient  recorder.  The  transient  recorder  accomplishes 
the  A/D  -  conversion  and  stores  the  digitized  time  record  in  its 
memory.  An  interface  bus  allows  the  computer  to  have  access  to  the 
seven  channels  (memories)  of  the  transient  recorder.  Further  data 
processing  is  done  under  the  control  of  the  program  listed  in  the 
Appendix. 

5.1.2.  This  program  is  written  in  BASIC  and  runs  on  the  desktop 
computer  Tektronix  4054.  The  three  parts  of  the  program  are  stored 
on  the  disc  in  three  files  (a  back  up  tape  cartridge  contains  the 
program  in  binary  code).  Part  I  reads  the  calibration  step  and  the 
time  record  from  the  memory  of  the  transient  recorder,  and  calculates 
shock  spectra.  Part  II  contains  all  output  procedures.  Part  III  is 
basically  the  same  as  Part  I  except  that  it  reads  a  calibrated 
digital  time  record  from  the  external  tape  unit  rather  than  from  the 
transient  recorder.  Both  Part  I  and  Part  III  use  the  output 
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procedures  in  Part  II. 

5.2  Starting  Up 

The  computer  and  at  least  half  of  the  units  connected  with  the 
computer  by  interface  bus  must  be  turned  on.  To  run  the  program,  set  the 
clock  on  the  disc  drive,  mount  the  unit,  load  the  program  and  enter  "RUN". 
Example:  Enter:  CALL  "SETTIM",  "20  -  MAR  -81  □  8:30" 

Press :  RETURN 

Enter:  CALL  "MOUNT", t  ,  A$ 

Press:  RETURN 

Enter:  OLD  "  @  SPRGM/PART  I" 

Press:  RETURN 
Enter:  RUN 
Press:  RETURN 

Message  on  the  Screen: 

SPRGM/PART  I  is  loaded  (default  conditions)! 

The  corresponding  procedure  for  using  binary  tape  is  either: 

Enter:  FIND  2 
Press:  RETURN 
Enter:  CALL  "BOLD" 

Press:  RETURN 
Enter:  RUN 
Press:  RETURN 
or  simply  press:  AUTOLOAD 

While  running  the  program,  the  computer  indicates  a  need  of  communication 
with  the  user  by  a  "beep"  tone  and  a  message  on  the  screen. 

The  symbol  Q  indicates  one  space. 

5.3  Default  Conditions 

Default  conditions  are  listed  in  line  290  to  500  in  the  Appendix.  All 
default  parameters  may  be  changed  at  any  time  by  entering  the  desired  values 
from  the  keyboard  (e.g.  entering  F2  *  500  changes  the  upper  analysis  frequency 
to  500  Hz).  The  current  value  of  every  parameter  can  be  checked  by  entering 
the  parameter  and  pressing  RETURN:  the  current  value  of  the  parameter  appears 
on  the  screen. 

The  most  important  parameters  defined  by  default  conditions  are  explained 
below : 

FI,  F2 :  Lower/upper  analysis  frequency:  can  be  any 

positive  value  between  1  and  1000. 

C0:  Number  of  analysis  frequencies  between  FI  and 

F2:  depends  on  memory  available  (Part  I  up  to 
about  110,  Part  III  about  160).  C0  »  1  is 
possible  also. 

F8 :  Cut  off  frequency  of  low  pass  filter  (default 

1000  Hz). 

Nl,  N2 :  Lower/upper  value  of  y-axis  in  log  display 
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N6 ,  N7  :  Lower/upper  value  of  x-axis  in  log  display 

D:  Primary  address  for  output  devices: 

32  screen 

3  ext*  tape 

4  plotter 
40  printer 

P:  selects  grid  in  log  display: 

0<P<10  sets  grid  line  at  P 
P>10  no  grid 
(P  must  not  be  zero!) 

0:  indicates  the  part  of  the  program  currently 

loaded  in  the  memory 

1.  Part  I 

2.  Part  II 

3.  Part  III 

(Do  not  change  parameter  0  from  the  keyboard) • 
M9:  Calibration  matrix. 

A2 :  Offset-correction  matrix. 

5.4.  User  Defined  Keys 


The  first  row  of  user  defined  keys  (No.  11  through  15)  are  used  with 
program  Part  I  or  III.  All  other  user  defined  keys  (No.  1  through  10  and  16 
through  20)  belong  to  program  Part  II.  User  Defined  Key  No.  4  and  _5_  allow 
switching  from  program  Part  I  or  III  to  Part  II  or  vice  versa.  Pressing  a 
user  defined  key  which  belongs  to  a  program  part  not  currently  loaded  in  the 
memory  causes  an  error  message  to  appear  on  the  screen. 

The  following  describes  each  user  defined  key  and  the  system's 

reaction. 


5.4.1.  First  Row  (No.  11  through  15) 

User  Defined  Key  No  11  triggers  the  calibration  routine  in  program 
Part  I.  After  the  time-record  has  been  defined  as  a  displacement, 
velocity  or  acceleration  record  there  are  two  ways  to  fill  the 
calibration  matrix  M9  with  the  calibration  factors  for  each  channel. 

First  the  calibration  factor  may  be  entered  from  the  keyboard  in  input 
units  per  1  Volt,  or  secondly,  the  calibration  step  may  be  read  from  each 
channel  of  the  transient  recorder.  In  the  latter  case  the  calibration 
level  has  to  be  entered  in  appropriate  units  before  the  computer  reads 
the  content  of  each  channel  and  displays  it  on  the  screen.  The  user  may 
define  the  calibration  level  by  means  of  setting  the  crosshair  first  to 
the  upper  and  then  to  the  lower  calibration  level.  The  calibration 
factor  (i.e.  input  unit/digit)  for  this  particular  channel  will  be 
assigned  to  an  element  of  calibration  matrix  M9.  The  program  starts  this 
procedure  all  over  again  until  the  calibration  steps  of  all  7  channels 
have  been  read.  If  the  program  encounters  a  channel  set  to  ground  or 
calibration  mode  an  appropriate  message  appears  on  the  screen  and  no 
calibration  factor  will  be  defined  for  this  particular  channel.  The 
program  proceeds  with  the  next  channel. 

User  Defined  Key  No.  12  reads  the  time  record  from  a  defined 
channel  of  the  transient  recorder  and  displays  a  plot  of  acceleration. 
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velocity  and  displacement's  time  on  the  screen.  Also  each  maximum  and 
minimum  value  is  printed.  When  running  program  Part  I,  the  calibration 
routine  (key  11)  must  be  completed  before  using  Key  No.  12.  With  program 
part  III,  User  Defined  Key  No.  12  causes  the  computer  to  read  the  time 
record  (including  other  important  information  like  calibration  factor, 
time  step,  etc.)  from  the  external  tape. 

The  read  routine  has  to  be  completed  for  each  record  before 
calculating  shock  data  by  pushing  User  Defined  Key  No.  13.  The  computer 
then  calculates  shock  data  of  the  time  record  that  has  been  input  by 
means  of  Key  1 2 . 

Using  the  mathematical  algorithm  described  above,  the  program 
computes : 

1.  Analysis  Frequency 

2.  Pseudo  velocity  Maximax-Shock  Spectrum 

3.  Rel.  Velocity  -  Residual- Shock  Spectrum  (or 
Fourier  Transform) 

4.  Rel.  Displacements 

5.  Abs.  Acceleration 

6.  Time  when  maximum  occurs 

7.  Phase 

While  calculating,  the  computer  displays  on  the  screen  the  serial 
number  of  analysis  frequency  (from  1  at  the  lowest  to  CJ&  at  the 
highest),  the  analysis  frequency,  the  pseudo-velocity ,  and  the 
number  of  data  points  used  for  the  calculation  of  shock  spectra  for 
each  analysis  frequency.  The  computing  time  is  also  given  at  the 
end  of  the  procedure. 

The  number  of  analysis  frequencies  (1  to  160  ),  as  well  as  any 
interesting  range  of  analysis  frequencies  between  1  and  1000  Fz,  can 
be  chosen.  (Note  that  F2  has  to  be  larger  than  Fj). 

The  time  record  and  all  important  information,  such  as 
calibration  factor,  time  step,  etc.  can  be  stored  on  a  defined  file 
of  the  external  tape  by  pushing  User  Defined  Key  No.  14.  Additional 
analysis  can  be  done  on  this  time  record  by  means  of  part  III  of  the 
program  (refer  to  User  Defined  Key  No.  12).  An  offset  correction  of 
the  time  record  can  be  made  with  the  procedure  of  User  Defined  Key 
No.  15.  The  first  integral  of  the  time  record  is  displayed  on  the 
screen  and  the  user  may  define  the  zero  line  by  means  of  the  cross 
hair.  After  pushing  any  key  the  computer  calculates  an  offset 
correction  for  this  particular  channel  and  considers  that  at  all 
later  operations.  This  procedure  may  be  repeated  several  times. 

*  The  title  or  description  of  the  time  record  is  stored  in  string  W$  and 
appears  on  the  top  of  the  display.  It  may  be  changed  anytime  by  redefining  W  $ 

^Depends  on  memory  available  for  the  calculation  (program  part  III  allows  up 
to  160  analysis  frequencies). 
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However,  the  user  has  to  realize  that  this  procedure  requires  some  actual 
value  to  be  known. 


5.4.2.  Second  Row  (No.  1  through  5) 

The  results  of  the  shock  spectra  computation  (Key  No.  13)  can 
be  saved  by  pushing  User  Defined  Key  No.  1.  This  causes  the 
computer  to  write  all  data  on  a  defined  file  of  the  external  tape. 
These  shock  data  can  be  read  back  with  User  Defined  Key  No.  2. 

A  numerical  output  of  all  analysis  data  can  be  retrieved  by 
pushing  User  Defined  Key  No.  3.  As  for  all  other  output  procedures 
the  primary  address  of  the  selected  output  device  (screen,  printer, 
or  plotter)  has  to  be  specified  here. 

User  Defined  Key  No.  4  makes  the  computer  load  part  II  of  the 
program  while  keeping  all  values  of  all  parameters  and  variables  in 
its  memory. 

5.4.3.  Third  Row  (No.  16  through  20) 

The  third  row  of  user  defined  keys  produces  a  linear  plot  of 
the  shock  sgectra.  The  amplitude  scale  is  automatically  adjusted  to 
the  output.  Any  range  of  analysis  frequencies  can  be  selected. 

User  Defined  Key  No.  16  displays  a  plot  of  pseudo-velocity 
vs.  analysis  frequency  (maximax) 

User  Defined  Key  No.  17  displays  a  plot  of  rel.  velocity 
vs.  analysis  frequency  (Residual) 

User  Defined  Key  No.  18  displays  a  plot  of  acceleration  vs. 
analysis  frequency . 

User  Defined  Key  No.  19  displays  a  plot  of  rel. 
displacements  vs.  analysis  frequency. 

User  Defined  Key  No.  20  displays  a  plot  of  phase  (0  to  180 
degree)  vs.  analysis  frequency. 

5.4.4.  Fourth  Row  (No.  6  through  10) 

A  logarithmic  plot  of  analysis  data  can  be  retrieved  by 
pushing  any  of  the  user  defined  keys  in  the  fourth  row.  User 
Defined  Key  No.  6  displays  the  Maximax  Shock  Spectrum,  No.  7  the 
Residual  Shock  Spectrum,  and  No .  8  both  of  them.  User  Defined  Key 
No.  9  allows  some  analysis  work  on  the  screen  in  determining  modal 
frequencies  and  appropriate  shock  spectrum  values  (the  Design  Shock 
Spectrum)  by  placing  the  two  crosshairs  at  points  of  interest.  The 
thumbwheels  on  the  keyboard  control  the  movement  of  the  crosshairs 
When  a  keyboard  key  is  pressed,  the  location  of  the  crosshair  is 
marked  and  the  appropriate  analysis  frequency  and  shock  spectrum 
value  is  printed  beside  the  graph-plot.  When  the  key  "0"  is 


* 


except  phase  is  always  displayed  between  0  and  180  degrees. 


pressed,  the  crosshair  disappears  from  the  screen.  (Note:  this 
procedure  works  with  the  screen  only!) 


User  Defined  Key  No.  10  gives  a  plot  of  a  number  of  different 
shock  spectra  stored  on  the  external  tape.  Each  spectrun  is  drawn 
with  a  different  dark  line  and  the  description  of  the  spectrum  is 
printed  wherever  the  crosshairs  have  been  placed. 

6.  Examples  of  Analysis  System  Output. 

6.1.  Figures  3  and  4  show  the  various  output  formats  which  can  be 
selected.  These  were  made  using  the  hard  copy  unit,  which  produces  a 
reduced-size  version  of  the  CRT  display.  Higher  quality  graphs  can  be 
obtained  from  the  plotter. 

6.2.  For  Fig.  3,  the  input  record  was  a  half-sine  acceleration  pulse. 

The  original,  integral  and  double  integral  are  shown  in  Fig.  3A  as  scaled 
acceleration,  velocity  and  displacement.  The  raaximax  shock  spectrum  is 
shown  in  linear  plot,  without  grid-lines,  in  Fig.  3B,  and  in  log-log 

^  plot,  with  grid-lines,  in  Fig.  3C*  The  raaximax  and  residual  are  shown  in 
log-log  plot,  with  factor-of-ten  grid-lines,  in  Fig.  3D.  Finally,  a 
tabular  listing  of  all  computed  data  is  shown  in  Fig.  3E. 

6.3.  The  input  record  for  Fig.  4  was  an  acceleration-t ime  history 
recorded  on  a  resilient  structure  aboard  the  Floating  Shock  Platform  with 
shot  standoff  of  20  ft.  The  acceleration,  velocity  and  displacement 
waveforms  are  shown  in  Fig.  4A,  utilizing  the  offset  correction  routine 
described  above:  i.e.,  requiring  the  average  velocity  near  the  400 
millisec.  mark  to  be  zero.  The  maxiraax  and  residual  spectra  are  given  in 
Fig.  4B,  and  again  in  Fig.  4C,  which  demonstrates  the  extraction  of 
design  shock  spectral  parameters.  Modal  frequencies  and  shock  spectrum 
values  are  lined  up  with  the  CRT  crosshairs,  and  their  values  tabulated. 
The  values  read  are  also  identified  on  the  graph.  The  tabular  listing  of 
computed  data  is  shown  in  Fig.  4D. 
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1 1««  record  (a)  Acceleration,  velocity  and  displacement  time  histories 

Fig.  3  —  Output  for  analysis  of  a  half-sine  acceleration  pulse 
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—  Output  for  analysis  of  a  half-sine  acceleration  pulse 
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Maximax  and  residual  shock  spectra,  log-plot,  with  factor-of-te 
Fig.  3  —  Output  for  analysis  of  a  half-sine  acceleration  pulse 
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(e)  Tabular  listing  of  all  computed  data 
Fig.  3  —  Output  for  analysis  of  a  half-sine  acceleration  pulse 
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(a)  Acceleration,  velocity,  and  displacement  time  histories 
Fig.  4  —  Output  for  analysis  of  FSP  acceleration  record 
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(b)  Maximax  and  residual  shock  spectra 
Fig.  4  —  Output  for  analysis  of  FSP  acceleration  record 
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(d)  Tabular  listing  of  all  computed  data 
Fig.  4  —  Output  for  analysis  of  FSP  acceleration  record 
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Appendix  1 

List  of  Program  Variables 


1  A  Array  for  rime  record  QIIKN0) 

2:A2  Zero  offset  -  Motrix  0111(7,2* 

3  0  Arroy  0  T M (2 ) 

4  03  Variable  used  in  subroutine  TIME 
504  Vorioole  used  in  suorour ina  TIME 
6  C  Arroy  OIM(C0) 

7 ; C0  Actual  number  of  analysis  frequencies  (default  90) 

8 : C 1  El emen t  of  F5  (first  ana  lysis  f requency  per  f r equency 
9 • C2  Element  of  F5  (Iasi  analysis  f requency  per  frequency 
10. C3  Max  possible  analysis  frequencies  (depends  on  memory 
<1:0  I/O  device  address  (default  32) 

12: 01  rile  number  on  ext.  tape  drive  to  store  data 

1303  File  number  on  ext.  tape  drive  to  read  data  from 

14.E  Instantaneous  factor  for  sample  step 

15  El  Factor  for  sample  step  within  frequency  range  1M-10HZ) 

16  E2  Factor  for  sample  step  within  frequency  range  IIU0-100Hz» 

17;E3  Factor  for  sample  step  within  frequency  ronge  1 1  I ( 1 00- 1 000Hz ) 

18;F  Instantaneous  value  of  analysis  frequency 

19F1  Lowest  analysis  frequency  (default  1) 

20:F2  Highest  analysis  frequency  (default  1000) 

21  ;F3  ^ote'*  analysis  frequency  within  one  of  the  3  frequency  ranges 

2 2:FA  Upper  analysis  frequency  within  one  of  the  3  frequency  ranges 

23:F5  Frequency  array  DIh(C0) 

24:F6  Lower  frequency  of  lin.  graph  (input) 

25=F7  Upper  frequency  of  lin.  graph  (input) 

26:F8  Cutoff  frequency  of  lowpass  filter 

27  G  Instantaneous  value  of  accelarat ion  ( (2*p i fF) f2*X) 

28  G5  Acceleration  array  0in(C3) 

29  H  Time  increment  of  digitized  time  record 
30: HI  Exponent 

3'  H2  Time  increment  for  computation 

32  I  Coup  ter 

33  It  Counter  to  find  Cl 

34  12  Indi cares  f requency  range  (I  to  III) 

3513  Sample  number  when  V2  occurs 

36.14  Sample  number  when  VI  occurs 
37  K  Target  or~oy  for  graph  disp  lay 

38'M  Factor  of  F8  to  calculate  sample  frequency  within  frequency  ronqe  I 
e* aul t  2)  * 

39  N2  Focror  of  F8  to  calculate  sample  frequency  within  frequency  range  II 
def cjIt  2) 

40  *3  Factor  of  F8  to  calculate  sample  frequency  wi th:n  frequency  range  II 
(def  aul t  1 0) 

41  L  Calibration  level  in  \ npu t  units 
42. M  Number  of  data  module 

43  Ml  Maximum  of  array  \  for  outoscole 

44  M9  Calibration  Matrix  0IM(7,2) 

45  N0  Number  of  data  points  on  record  (default  4096) 

46N  Njmber  of  data  points  considered  dur  i  ng  calculator 

47  N1  _ower  vert >001  level  of  log  graph  display 

48  N2  Upper  vert. col  level  of  log  gropn  display 

49. N4  Number  of  data  points  at  the  beg  i  nn  i  ng  of  rne  residual  i-esDorse 
50  •  N5  Number  of  data  points  after  fruacror ,on 
5’  N6  Lower  frequency  of  log  graph  display 
52  N?  Upper  f requency  of  log  graph  display 

53. C  Indicates  *he  par t  of  the  program  currently  loaded  M  or  21 
54  P  Selects  g~ - d  on  log,  display  fP>!0  no  grid,  default  10) 


r angel 
ronge) 
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55. PI  Error  indicator  if  >0  902 

56  Q  Number  of  periods  considered  oer  analysis  frequency 

57 . R  Selects  MAX  I  MAX-  and/or  RESIDUAL  SHOCK  SPECTRun ( Keys  6  to  9l 

58  R1  Selects  d»v.  linear  spectra  (user  defined  Keys  16  to  22) 

59  R2  Ind  cates  di splacemenr ( 1 ) , veloci ry (2)  or  accel erot i on ( 3 )  record 

60  R3  Equals  X3  before  N4  (initial) 

6*.-R4  Equals  X3  after  N4  (residual! 

62  R5  Flag  for  of f se t -procedure 

63  R6  Flag  for  deleting  calibration  rout  ire 
©4  R8  Flag  in  calibration  routine 

65  S  Calibration  factor 

66  S 1  Calibration  factor  after  di f  f  erenr iot ion  and/or  integration 
67:S2  Calbrotion  factor  after  E-procedure 

63  T  I  ns t on t aneous  value  of  time  of  maximum  ( T  3  or  T4) 

69  T t  Beg i nn i ng  of  time  record 

70  T3  Time  of  max.  response  before  N4  (inifiol) 

7 '  t 4  Time  o£  max.  response  after  N4  (residual) 

72  T5  Time  array  (contO'ns  max  values  of  either  T3  or  T4)  OIM(C0) 

73  U  Time  difference  m  decimal  value 

7 4  V  Instonroneous  value  of  pseudve 1 oc i t y  ((2*pi*F)*X) 

75 : V •  Max i mum  of  orray  A 

76  V2  Mini  mum  of  array  A 
77 : v5  Veloci»y  array  OIH<C0) 

78 ; w  Foe  t or 

7g  y •  Number  of  spectra  pro  plot  (Multi  Display  Mode)  and  flag  for  wt-qispio 

y 

90  V2  Array .  DIM  W2(W1 ) 

81  X  X-coordmate  of  pointer 

82  X  Instantaneous  value  of  displacement  (maximum  of  R3  or  R4 i 
63  X-  Deflection  at  -0.5*h2 

84 : X2  Oeflection  at  +0.5SH2 

65;  X3  Max.  response  of  analysis  frequency  F 

86:X5  Displacement  array 

87 ; v  v-coordinote  of  pointer 

88  Y  • ns r on r aneous  value  of  of  mox .  response  velocity  after  N4 
8S: Y5  Velocity  orroy  OIM(C0) 

90 ; Z  I rs t ant aneous  value  of  of  phase 
9i:Z5  Phase  orray  Drfl(C0) 

92;  AS  String  for  graph  display 

93 = 3*  System  time  when  process-procedure  is  completed 
94 : CS  Time  string 

95 . DS  String  :o  address  data  module 
96  GS  String  for  g^oph  display 
97 ; HS  String  for  graph  display 

98  Li  String  *0  checK  cond. ’ ion  of  data  module 

99  MS  String  for  graph  display 
J 00  NS  String  for  graph  display 

1 0 1  PS  Spring  for  graph  display 

102  OS  Spring  (  f  *0"  stops  analysis  procedure) 

03:RS  "NPL  -  CODE  5842" 

104  TS  Contains  system  time 
05 : V s  Contains  text  for  vert  level 
’  06  VS  Title  to  identify  record 
I  07; XS  Segment  spring 
'08  xs  Segmen t  string 
109  Zt  Segment  string 


Appendix  2 
Program  Listing 

1  REM  ////////////////////////  8SPRGM/PART  I  \\\\\\\\\\\\\\\\\\\\\\\\\\\ 

2  INIT 

3  GO  TO  260 

4  IF  0=1  THEN  92 

5  GO  TO  7450 

8  IF  0=1  THEN  92 

9  GO  TO  7530 

12  IF  0=1  THEN  92 

13  GO  TO  8980 
16  GO  TO  190 
20  GO  TO  100 

24  IF  0=1  THEN  92 

25  R*1 

26  GO  TO  7690 

28  IF  0=1  THEN  92 

29  R=2 

30  GO  TO  7690 

32  IF  0=1  THEN  92 

33  R  =3 

34  GO  TO  7690 

36  IF  0=1  THEN  92 

37  GO  TO  9510 

40  IF  0=1  THEN  92 

41  GO  TO  9780 

44  IF  0=2  THEN  90 

45  CHARSIZE  3 

46  IF  R6=0  THEN  94 

47  IF  R6= 1  THEN  4200 

48  IF  0=2  THEN  90 

49  GO  TO  530 

50  DELETE  A 

52  IF  0=2  THEN  90 

53  GO  TO  2510 

54  DELETE  A 

56  IF  0=2  THEN  90 

57  GO  TO  5250 

60  IF  0=2  THEN  90 

61  R5= 1 

62  CHARSIZE  3 

63  GO  TO  2410 

64  IF  0=1  THEN  92 

65  R 1 =1 

66  GO  TO  6010 

68  IF  0=1  THEN  92 

69  R 1 =4 

70  GO  TO  6010 

72  IF  0=1  THEN  92 

73  R 1 =2 

74  GO  TO  6010 

76  IF  0=1  THEN  92 

77  R 1 =3 

78  GO  TO  6010 

80  IF  0=1  THEN  92 

81  R 1 =5 

82  GO  TO  6010 

90  PRINT  "You  qot  to  PRESS  KEY  5  first 'G" 

91  RETURN 

92  PRINT  "You  qot  to  PRESS  KEY  4  first'G" 

93  RETURN 

94  PRI  "LJJTHE  CALIBRATION  PROCEDURE  HAS  BEEN  DELETED  IN  ORDER  TO  SAUE"; 
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95  PRINT  M  MEMORY.  IF  YOU  WANT  TO  REDEFINE  THE  CALIBRATION  MATRIX  M9  "  ; 

96  PRINT  "RELOAD  THE  PROGRAM  (PRESS  USER  DEFINED  KEY  4  AND  5  IN  A  ROW" 

97  PRINT  "OR  PRESS  > AUTOLOAD< ) * GG  " 

99  RETURN 

100  REM  ***************  LOAD  ROUTINE  *********************************** 
110  IF  0= 1  THEN  160 

120  DELETE  F5 , X5 , VS , G5 , T5 , Y5 . 25 , C , K *  A 
130  R 6= 1 
140  0=1 

150  APPEND  "@SPRGM/P ARTI " ; 1 

160  PRINT  "LGSPRGM/PART  I  IS  LOADED!" 

170  DELETE  5380,20000 
130  RETURN 
190  IF  0=2  THEN  240 
200  0=2 

210  DELETE  270,20000 

220  APPEND  "8SPRGM/ PARTI I " ; 260 

230  RENUMBER  6000,10,260 

240  PRINT  "LGSPRGM/PART  II  IS  LOADED!" 

250  RETURN 

260  REM  ***************  DEFAULT  CONDITION  ****************************** 
270  W$  =  "TEST " 

280  0=1 
290  W  1  =  0 
300  R5=0 
310  R6= l 
320  D=32 
330  F 1 = 1 
340  F2= 1 000 
350  F8= 1000 
360  N 1 =  0 , 0 1 
370  N2= 1 0  0 
380  N6=l 
390  N7= 1000 
400  CQ=90 
410  N0=4096 
420  K 1=2 
430  K  2  =  2 
440  K  3  = 1 0 
450  P=10 

460  DIM  M9(7,2) ,C(C0 ) ,A2<7,2> ,R$< 15) ,T*( 18) ,B$( 18) ,C$( 18) 

47  0  DIM  D*<8> ,Ht<20>  »Vt<22) ,Gt<7> ,M*<7) >Nt<8) ,X«<22> ,P*<35> ,Lt<4> ,Qt< 1 > 
480  R*="NRL  -  CODE  5842" 

490  M 9 r  0 
500  A2=  0 

510  PRINT  "  GSPRGM/P ART  I  IS  LOADED  (DEFAULT  CONDITIONS)!" 

520  RETURN 

530  REM  ******************  READ  RECORD  ********************************* 
540  DELETE  4200,5240 
550  R 6  =  0 

560  PRINT  @D,17:1.2,1.3 
570  CHARSI2E  3 

580  PRI  "LTHE  CURRENT  TITLE  IS  J  ENTER  A  NEW  TITLE  IF  NECESSA" ; 

590  PRINT  "RY.  OTHERWISE  PRESS  'RETURN '  . JG " 

600  INPUT  A$ 

610  IF  A*="M  THEN  630 
620  W4=AS 
630  DIM  A  <  NO ) 

640  CHARSI2E  3 

650  IF  R5=  0  THEN  660 


28 


660  W3* 1 30 

670  IF  D< >4  THEN  690 
680  U3«150 

690  PRINT  “ JENTER  NUMBER  OF  DATA  MODULE  TO  READ  FROM ;  G  “ ; 

700  INPUT  M 

710  IF  R5=*l  THEN  740 

720  IF  M9<  M , 2 ) =0  THEN  2390 

730  S=M9 ( M ,  2 ) 

740  PRINT  96: "P" 

750  GO  TO  M  OF  760,780,800,820,840,860,880 
760  D$a "101 B4096 " 

770  GO  TO  890 
780  D*= " I 02B4096 " 

790  GO  TO  890 
800  D$= " I 03B4096 " 

810  GO  TO  890 
820  D*= ”10 4B4096 " 

830  GO  TO  890 
840  D$= " I 05B4096 " 

850  GO  TO  890 
860  D*= " I 06B4096 “ 

870  GO  TO  890 
880  D$=  M 1 07B4096 " 

890  PRINT  86:D$ 

900  INPUT  86 : A 
910  P*="DM 
920  PRINT  86  :  P$ 

930  INPUT  86  :  P* 

940  A=*A-A2  ( M  ,  2 ) 

950  H$=“SA“ 

960  PRINT  86:H* 

970  INPUT  86 : H$ 

980  H1=VAL <H*> 

990  IF  HI >12  THEN  1180* 

1000  GO  TO  HI  OF  1 , 1 , 1 , 1 0 1 0 , 1 030 , 1 050 , 1 070 , 1 090 , 1 1 1 0 , 11 30 , 11 50 , 1 1 70 

1010  H=*0  .  005 

1020  GO  TO  1200 

1030  H=0  ,  0 1 

1040  GO  TO  1200 

1050  H=0.025 

1060  GO  TO  1200 

1070  H=  0  .  05 

1080  GO  TO  1200 

1090  H=0 . 1 

1100  GO  TO  1200 

1110  H=0 .25 

1120  GO  TO  1200 

1130  H*0 . 5 

1140  GO  TO  1200 

1150  H* 1 

1160  GO  TO  1200 

1170  Ha*2 . 5 

1180  PRINT  “ENTER  TIME  INCREMENT  IN  MS 

1190  INPUT  H 

1200  IF  R5=l  THEN  2430 

1210  CHARSIZE  2 

1220  CALL  “TIME  “  , T* 

1230  PRINT  8D; ”LK" jR4 
1240  PRINT  8D:T$ 

1250  WINDOW  0,W3, 0,100 
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1260  VIEWPORT  0>U3, 0,100 

1270  IF  W1 >0  THEN  1350 

1280  PRINT  8D, 17 : 2,3 

1290  CHARSIZE  4 

1300  MOVE  BD : W3/2 , 98 

1310  FOR  1=1  TO  LEN(W*>/2 

1320  PRINT  BD: "H" 1 

1330  NEXT  I 

1340  PRINT  BD : Wt 

1350  PRINT  BD , 1 7 : 1 .2,1 .8 

1360  CHARSIZE  2 

1370  T1=H*ABS(VAL(P*) > 

1380  GO  TO  R2  OF  1390,1520,1660 

1390  GOSUB  1840 

1400  S1=S 

1410  GOSUB  2020 

1420  PRINT  BD : “ * " 

1430  CALL  "DIF3 11 , A , A 

1440  SI =51 /H 

1450  GOSUB  1900 

1460  GOSUB  2020 

1470  CALL  "DIF3" , A, A 

1480  S1=S1 *10A3/(H*9.8044> 

1490  GOSUB  1960 
1500  GOSUB  2020 
1510  GO  TO  1780 
1520  GOSUB  1900 
1530  S1=S 
154 0  GOSUB  2020 
1550  PRINT  BD:"*" 

1560  CALL  "DIF3",A,A 
1570  SI =S 1*10*3/(9.80 44*H> 

1580  GOSUB  1960 

1590  GOSUB  2020 

1600  CALL  "INT  ”,A,A 

1610  CALL  "INT  ",A,A 

1620  S1=S1*1 . 0E-3*HA2»9.8044 

1630  GOSUB  1840 

1640  GOSUB  2020 

1650  GO  TO  1780 

1660  GOSUB  1960 

1670  S1=S 

1680  GOSUB  2020 

1690  PRINT  BD;"*" 

1700  CALL  "INT", A, A 

1710  S1=S1*1 .0E-3*H»9.8044 

1720  GOSUB  1900 

1730  GOSUB  2020 

1740  CALL  "INT", A, A 

1750  S1=S1*H 

1760  GOSUB  1840 

1770  GOSUB  2020 

1780  VIEWPORT  0,W3, 0,100 

1790  WINDOW  0,W3, 0,100 

1800  MOVE  BD : 0 , 0 

1810  PRINT  BD:"*  oriq.  tine  record 
1820  DELETE  A 
1830  RETURN 

1840  VIEWPORT  5,W3-30,4,32 
1850  VS* "DISPLACEMENT" 


I860  G*=“D  nax  j  M 
1870  M*="D  Min  :  " 

1880  N$="  [mm3" 

1390  RETURN 

1900  VIEWPORT  5,W3-30,36,64 
1910  V$  = "VELOCITY " 

1920  G*="V  nax  : " 

1930  M$="V  win  : " 

1940  N$="  [d/sec]" 

1950  RETURN 

1960  VIEWPORT  5 , W3-30 , 68 , 96 
1970  V$= "ACCELERATION " 

1980  G$="A  Max  : " 

1990  H$*"A  Min  : " 

2000  N$= "  C q ] " 

2010  RETURN 

2020  CALL  "MAX", A, VI ,14 
2030  CALL  "MIN " , A , V2 , 13 
2040  WINDOW  0 , NO , V2 , VI 
2050  IF  D=32  THEN  2110 
2060  MOVE  @D : 0 , 0 
2070  FOR  1=1  TO  NO 
2080  DRAW  8D : I >  A  < I ) 

2090  NEXT  I 

2100  GO  TO  2120 

2110  CALL  "D ISP  ",  A 

2120  AXIS  0D:4OO,O 

2130  IF  R5=l  THEN  2340 

2140  FOR  1=0  TO  NO  STEP  2000 

2150  MOVE  @D: 1,0 

2160  PRINT  @D :  "HJ" ,l*H\ 

2170  NEXT  I 

2180  PRINT  3D: "HHHHHHHHHHH C MS 3 " 

2190  MOVE  @D : 4200 , VI 

2200  PRINT  @D :  USING  2240 : " J " ,G$ , VI «S1 , N* ; 

2210  MOVE  @D : 420  0  , VI 
2220  PRINT  @D:  '*  J  r'  j 

2230  PRINT  @D :  USING  2240 i " J M , M* , V2*S1 , N$ > 

2240  IMAGE  A  *  7A  >  6D , 1  D  >  1 1 A 
2250  MOVE  80 ; 0 , 0  *  5*  < V2+V1 ) 

2260  PRINT  8D:  "HH"  ; 

2270  FOR  1=1  TO  LEN<V*>/2 
2280  PRINT  @D :  "K"; 

2290  NEXT  I 

2300  FOR  1=1  TO  LEN<V*> 

2310  X*=SEG<V$, I  ,  1  ) 

2320  PRINT  @D;X*; "HJ"; 

2330  NEXT  I 
2340  RETURN 

2350  PRINT  @D:"JDATA  MODULE  "jMj"  IS  IN  GROUND  MODE  MGG” 
2360  RETURN 

2370  PRINT  " JDATA  MODULE  "  ;  h ;  "  IS  IN  CALIBRATION  MODE  ! ! GG " 
2380  RETURN 

2390  PRINT  " JREAD  THE  CALIBRATION  STEP  FIRSTIG" 

2400  RETURN 

2410  PRINT  "LQFFSET  -  CORRECTION:" 

2420  GO  TO  630 
2430  VIEWPORT  0,110,10,90 
2440  CALL  "INT",A,A 
2450  GOSUB  2020 
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2460  POINTER  X,Y,Q* 

2470  A2<M,1)=M 

2480  A2  <  M , 2 ) =Y/X+A2 <  M  >  2  > 

2490  R5=0 
2500  RETURN 

2510  REM  a*******************  PROCESS  *****#****x*********************** 
2520  CHARSIZE  3 
2530  CALL  "TIME " , B* 

2540  PRINT  @D :  "L " ; B$  j  "  J 11 
2550  PRINT  @32,26:2 
2560  IF  F 1 < 1  THEN  2620 
2570  IF  FI >1000  THEN  2620 
2580  IF  F2< 1  THEN  2630 
2590  IF  F2> 1 000  THEN  2630 
2600  IF  F2<  =F1  THEN  2630 

2619  GO  TO  2660 

2620  IF  C0  =  1  THEN  2660 

2630  PRI  “ANALYSIS  FREQUENCY  IS  OUT  OF  RANGE  ENTER  NEW  VALUES  OF  F1,F2G" 
2640  INPUT  F1,F2 
2650  GO  TO  2530 

2660  DELETE  F5 , X5 , V5 „ G5 , T5 , Y5 , Z5 , C > K , A 
2670  DIM  A ( NO ) 

2680  C3=INT< ( MEMORY/8-1 8 ) /8 > 

2690  DIM  F5<C0) , X5<C0>,V5<C0>  , G5  <  CO ) , T5(CQ>,Y5(C0) , Z5 ( CQ ) , C  <  CO ) 

2700  E1»INT<1Q00/<K1*F8*H)) 

2710  E2*INT (  1  0  00/ < K2*F8*H ) ) 

2720  E3=INT (1000/(K3*F8*H) > 

2730  F5 < 1 >=F 1 

2740  IF  C0=1  THEN  2780 

2750  FOR  1=2  TO  CO 

2760  F5( I )=<F2/F1 >A(1/(C0-1 ) )*F5<I-1 > 

2770  NEXT  I 
2780  X5=Q , 0 1 
2790  V5=0 . 0 1 
2800  G5=*0  .  0 1 
2810  T5=Q ,01 
2820  Y5=0 .01 
2830  Z5=0.01 
2840  11*1 
2850  12=0 

2860  IF  F 1 < 1 0  THEN  2900 
2870  12=1 

2880  IF  F 1 < 1 0 0  THEN  2900 

2890  12=2 

2900  GO  TO  3070 

2910  E=E1 

2920  G0SU6  3160 

2930  IF  II >C0  THEN  3010 

2940  GO  TO  3070 

2950  E*E2 

2960  GOSUB  3160 

2970  IF  I1>C0  THEN  3010 

2980  GO  TO  3070 

2990  E*E3 

3000  GOSUB  3160 

3010  DELETE  A 

3020  CALL  "TIME"  ,T$ 

3030  GOSUB  4030 
3040  PRINT  @32,26:0 

3050  PRINT  #D : "  MAX,  POSSIBLE  ANALYSIS  FREQUENCIES:  "  ;C3 


3060  RETURN 
3070  F3=F5< II ) 

3080  C1=I1 

3090  12=12+1 

3100  FOR  11=1  TO  CO 

3110  IF  F5  <  1 1 ) >10AI2  THEN  3130 

3120  NEXT  II 

3130  F4=F5(I1-1) 

3140  GO  TO  12  OF  291 0 > 2950 , 2990 

3150  RETURN 

3160  PRINT  @6 ; HPMjD$ 

3170  INPUT  86: A 
3180  A=A-A2(M>2) 

3190  IF  E<  =  1  THEN  3260 
3200  J=1 

3210  FOR  1=1  TO  NO  STEP  E 
3220  A  <  J ) =A  < I > 

3230  J=J+1 
3240  NEXT  I 
3250  GO  TO  3270 
3260  E=1 
3270  N5=N0/E 
3280  H2=H*E 

3290  GO  TO  R2  OF  330 0 > 3340 , 3370 

3300  CALL  " D I F 3 " >  A  >  A 

3310  CALL  "DIF3"  >A,A 

3320  S2=S 

3330  CO  TO  3380 

3340  CALL  "DIF3" ,A,A 

3350  S2=S*H2 

3360  GO  TO  3380 

3370  S2=S*1 0 A -3*H2A2*9 .8044 

3300  C2=I 1-1 

3390  Pl=0 . 0062831 853*F3*H2 
3400  FOR  I=C1  TO  C2 
3410  C  < I ) =2*C0S  <  P 1 ) 

3420  IF  Cl =C2  THEN  3440 

3430  P1=P1*<F4/F3)A(1/<C2-C1 )  ) 

3440  NEXT  I 

3450  IF  P 1 <  0 . 902  THEN  3500 
3460  IF  1 43 , 566/H2 >F2  THEN  3500 

3470  PRINT  "LGWARNING ;  FREQUENCIES  ARE  TOO  HIGH  FOR  YOUR  TINE  INCREMENT " 
3480  PRINT  ’*  ERRORS  ARE  OUTSIDE  THE  RANGE  -7  TO  *3  PERCENT  FOR " 

3490  PRINT  -  FREQUENCIES  ABOv>E  ;  1  43 . 566/H2  ;  "  HERTZ" 

3500  FOR  I=C1  TO  C2 
3510  F=F5(I) 

3520  IF  F<=1  THEN  3550 
3530  N=1000/<H2*5*LGT(F) ) 

3540  IF  N< N5  THEN  3560 

3550  N=N5 

3560  X1=0 

3570  X2=0 

3580  R3*0 

3590  R4-0 

3600  IF  F*0  THEN  3630 
3610  N4=N- 1-500/ < F*H2  > 

3620  IF  N4 >3  THEN  3640 
3630  N4*3 

3640  FOR  J=1  TO  N-l 
3650  X3=C(I)*X2-X1+A< J) 
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3660  IF  J<N4  THEN  3710 

3670  IF  ABS ( X3) < ABS < R 4 )  THEN  3740 

3680  R4  =  X3 

3690  T 4= J 

3700  GO  TO  3740 

3710  IF  ABS <  X3  >  <  ABS ( R3 )  THEN  3740 

3720  R3  =  X3 

3730  T3= J 

3740  XI =X2 

3750  X2=X3 

3760  NEXT  J 

3770  IF  ABS ( R3 ) > ABS<R4 )  THEN  3800 

3780  R3=R 4 

3790  T3=T 4 

3800  X-S2*R3 

3810  V=0 , 00628318*F*X 

3820  G-6 . 28831 8*F*V/9 .8044 

3830  T-Tl+(T3-0 , 5 ) *H2 

3840  Y=G . 00628318*F*S2*R4 

3850  Z=1 . 0E-3*<Tl+<T4-0 ,5)*H2)*F+0 ,25 

3860  Z=Z~INT(Z) 

3870  IF  Y  >  0  THEN  3890 
3880  Z-Z-0.5 

3890  IF  Z<  =0 . 5  THEN  3910 

3900  Z=Z-1 

3910  Z=360*Z 

3920  F5 < I ) =F 

3930  X5 ( I ) =X 

3940  V5< I >=V 

3950  G5  < I ) =G 

3960  T5 < I ) =T 

3970  Y5( I >=Y 

3980  Z5 ( I ) =Z 

3990  PRINT  @D:  USING  40 0 0 : I , F , V , N ; 

4000  IMAGE  20T,3D,4X,4D. 1 D ,4X , 5D . 2D , 3X , 6D 
4010  NEXT  I 
4020  RETURN 

4030  REM  ******  SUBROUTINE  TIME  ************************************** 
4040  DIM  B<2) 

4050 

4060  FOR  1*1  TO  2 
4070  X**SEG<C$, 1 1 ,8) 

4080  r$*SEG<XS>4,2) 

4090  Z$=SEG(X$,7 ,2) 

410  0  B< I >=VAL(X*>  +  <VAL( Y*)*60+VAL(Z*>  >/360  0 
4110  C$=T$ 

4120  NEXT  I 
4130  U*B ( 2 ) -B ( 1 ) 

4140  B3*(U~INT ( U ) )*60 
4150  B4=(B3-INT<B3>)*60 

4160  PRINT  @D  : * _ CALCULATION  TIME  (hr : Min : sec ) ;  •; 

4170  PRINT  8D:  USING  4 1 80  :  INT  ( U  >  ,  11 :  '\INT  ( B3 )  ,  H  ;  w  ,  INT  ( B4 ) 

4180  IMAGE  2D , A , 2D , A , 2D 
4190  RETURN 

4200  REM  ***************  CALIBRATION  *********************************** 
4210  DIM  A  ( NO  > 

4220  CHARSIZE  3 

4230  PRINT  "LARE  YOU  DEALING  WITH  J  A  DISPLACEMENT  RECORD  (MM)"; 

4240  PRINT  "  (ENTER  1)  A  VELOCITY  RECORD  (M/SEC)  (ENTER  2)  OR  A”  *> 

4250  PRINT  "  ACCELERATION  RECORD  (G)  (ENTER  3)g  -5 


4260  INPUT  R2 

4270  PRI  " _ DO  YOU  WANT _  TO  READ  A  CALIBRATION  STEP  (ENTER  1)  OR"; 

4280  PRINT  M  TO  ENTER  CALIBRATION  FACTOR  (ENTER  2>G  " ; 

4290  INPUT  R8 
4300  M9=0 

4310  IF  R8=2  THEN  4960 
4320  FOR  H- 1  TO  7 

4330  PRI  "LENTER  CALIBRATION  LEVEL  FOR  CHANNEL  " \H\"  IN  INPUT  UNITSG:" ; 

4340  INPUT  L 

4350  GOSUB  4750 

4360  IF  S4=0  THEN  4390 

4370  IF  S4=7  THEN  4420 

4380  GO  TO  4450 

4390  GOSUB  2310 

4400  CALL  "WAIT" ,2 

4410  GO  TO  4730 

4420  GOSUB  2330 

4430  CALL  “WAIT " , 2 

4440  GO  TO  4730 

4450  PRINT  @6: "P" 

4460  GO  TO  h  OF  447  0 ,4490  >4510  >4530  >4550 , 457  0  >4590 
4470  Df=  "10164096" 

4480  GO  TO  4600 
4490  D$= " I 02B4096 " 

4500  GO  TO  4600 
4510  D$= " I 03B4096 " 

4520  GO  TO  4600 
4530  D$= " I 04B4096  " 

4540  GO  TO  4600 
4550  D$= " IQ5B4096" 

4560  GO  TO  4600 
4570  D$=  M I 06B4096 " 

4580  GO  TO  4600 
4590  " 1 07B4096 " 

4600  PRINT  @6:D$ 

4610  INPUT  @6 : A 

4620  WINDOW  0,5000,-512,512 

4630  VIEWPORT  0,130,0,100 

4640  PRINT  "LCHANNEL  " ; M j  "  JJSET  CROSSHAIR  TO  UPPER  CALIBRATION  LEVELG M 
4650  CALL  "DISP",A 
4660  AXIS 

4670  POINTER  XI  ,Yt  ,Q$ 

4680  PRINT  "AJJJ  AND  LOWER  CALIBRATION  LEVELG!" 

4690  POINTER  X27Y2,Q* 

4700  IF  Y2 > Y 1  THEN  5230 
4710  H9  <  M , 1 ) 

4720  M9 (M,2)*L/(Y1-Y2) 

4730  NEXT  M 
4740  RETURN 

4750  GO  TO  M  OF  4760,4780,4800,4820,4840,4860,4880 
4760  L$= " A0 1 " 

4770  GO  TO  4890 
4780  L  $  =  "  A  0  2  " 

4790  GO  TO  4890 
4800  Li  =  " A03 " 

4810  GO  TO  4890 
4820  L $ - " A 0 4 " 

4830  GO  TO  4890 
4840  Lt="A05" 

4850  GO  TO  4890 


35 


4860  L$  =  "  A06 " 

4870  GO  TO  4890 
4880  L$="A07“ 

4890  PRINT  @6 ; Li 
4900  INPUT  @6:L$ 

4910  St=INT(VAL<L$)/l 000) 

492  0  S2= I NT < < UAL < L$ > ~S1 * 1 0 0 0 > / 1 0 0 > 

4930  S3=INT( (UAL(LS)-Sl *1  0QQ-S2*1 00 VI 0  ) 

4940  S4=VAL(L*)-S1*1000-S2*100-S3*10 

4950  RETURN 

4960  PRINT  "LENTER;" 

4970  FOR  M=1  TO  7 

4980  GOSUB  4750 

4990  IF  S4=0  THEN  5020 

5000  IF  S4=7  THEN  5040 

5010  GO  TO  5060 

5020  GOSUB  2310 

5030  GO  TO  5200 

5040  GOSUB  2330 

5050  GO  TO  5200 

5060  PRI  " JCAL IBR AT  ION  FACTOR  ( INPUT-UNITS/ 1 UGLT >  FOR  CHANNEL  M ; M j '* :  ” > 

5070  INPUT  L 

5080  GO  TO  S4  OF  5090 , 5 11 0 , 5 1 30 > 5 1 50 , 5 1 70 , 51 90 

5090  M9<M,2)=L/512 

5100  GO  TO  5200 

5110  M9 ( M , 2 ) =L*0 .2/512 

5120  GO  TO  5200 

5130  M9(M,2)=L *5/512 

5140  GO  TO  5200 

5150  M9(M,2)=L*0 . 1/512 

5160  GO  TO  5200 

5170  M  9  <  M  ,2)=L* 2/512 

5180  GO  TO  5200 

5190  M9(M,2)=L*0 .5/512 

5200  M9<M,1 )=M 

5210  NEXT  h 

5220  RETURN 

5230  PRINT  ,,AJ JJJJ  UPPER  CALIBRATION  LEUEL  FIRST!!  TRY  IT  AGAIN!  !G" 

5240  GO  TO  4670 

5250  REM  ***********************  WRITE  RECORD  ************************** 
5260  DIM  A ( NO ) 

5270  CHARSIZE  3 

5280  PRINT  "LENTER  FILE  NUMBER  (DEUICE  NR, 3)  TO  WRITE  RECORD  ONG s 
5290  INPUT  D3 
5300  FIND  @3:D3 
5310  PRINT  @6: 

5320  INPUT  @6  :  A 
5330  A=A-A2 (  M  ,  2  ) 

5340  WRITE  03:W$,R2,H,T1  >  M , M9 ( M , 2 ) , A2 ( M  ,  2 ) >  A 
5350  PRINT  " JRECORD  IS  STORED  IN  FILE  #G " ; D3 
5360  DELETE  A 
5370  RETURN 


6000  REM  SPRGM/P AR T  II  FILE  3 

6010  REM  ***********************  LIN.  SPECTRA  ************************** 
6020  DELETE  K 
6030  DIM  K ( CO  ) 

6040  AS= " MAXI MAX  SHOCK  SPECTRUM" 

6050  GO  TO  R 1  OF  6060,6080,6100,6120,6150 
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ram* 


6060  K=ABS ( V5) 

6070  GO  TO  6170 
6080  K=*ABS<  G5 ) 

6090  GO  TO  6170 
6100  K*ABS<  X5 ) 

6110  GO  TO  6170 

6120  A4= "RESIDUAL  SHOCK  SPECTRUM” 

6130  K=ABS(Y5> 

6140  GO  TO  6170 
6150  A4* “PHASE ” 

6160  K=Z5 

6170  CALL  ”MAX“,K,M1 ,14 
6180  CHARSIZE  3 

6190  PRINT  ML JFREQUENCY  RANGEG :  ” ; 

6200  INPUT  F6,F7 
6210  IF  D=4  THEN  6240 

6220  VIEWPORT  8*1.79,130-4*1.79,3*2.82,100-3*2.82 
6230  GO  TO  6260 

6240  VIEWPORT  8*1.79,150-4*1.79,3*2.82,100-3*2.82 
6250  PRINT  @D , 1 7 : 1 .2,1 .8 
6260  CALL  “TIME" ,T* 

6270  CHARSIZE  1 
6280  PRINT  8D: "L";R$ 

6290  PRINT  @D;T$ 

6300  IF  R 1 =5  THEN  6510 

6310  IF  ABS(M1 ><10  THEN  6350 

6320  IF  ABS ( Ml ) < 1 0 0  THEN  6390 

6330  IF  ABS ( hi ) <  1 00 0  THEN  6430 

6340  IF  ABS (Ml >>1000  THEN  6470 

6350  M1=«INT<M1  >  +  l 

6360  WINDOW  F6,F7,0,M1 

6370  AXIS  8D : <  F7-F6 ) /I 0,0.5 

6380  GO  TO  6580 

6390  M1=INT(M1/10)*10+10 

6400  WINDOW  F6,F7,0,M1 

6410  AXIS  @D : ( F7-F6 ) /I  0 , 5 

6420  GO  TO  6580 

6430  M1=INT(M1/100>*100+100 

6440  WINDOW  F6,F7,0,M1 

6450  AXIS  @D: ( F7-F6 ) /I 0 , 50 

6460  GO  TO  6580 

6470  Ml-I NT (Ml/1000)* 1000+1000 

6480  WINDOW  F6,F7,0,M1 

6490  AXIS  @D: ( F7-F6 ) / 1 0 , 50 0 

6500  GO  TO  6580 

6510  WINDOW  F6,F7, 0,180 

6520  AXIS  @D: < F7-F6 ) / 1 0 , 45 

6530  MOVE  £D ; F6 , 180 

6540  DRAW  @D:F7,1S0 

6550  DRAW  ®D:F7,G 

6560  MOVE  @D : (F7+F6)/2,180 

6570  CO  TO  6610 

6580  MOVE  0D : F6 , Ml 

6590  DRAW  fD:F7,Ml 

6600  DRAW  9D;F7,0 

6610  PRINT  8D, 17:2,3 

6620  CHARSIZE  4 

6630  GOSUB  6700 

6640  IF  W1 >0  THEN  6800 

6650  A*«U* 
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6660  PRINT  @D , 1 7 : 1  ,2,1 .8 

6670  CHARSIZE  2 

6680  GOSUB  6700 

6690  GO  TO  6800 

6700  IF  R 1 =5  THEN  6730 

6710  MOVE  &D: <F7+F6)/2,M1 

6720  GO  TO  6740 

6730  MOVE  @D: <F7*F6)/2, 180 

6740  PRINT  @D  :  “KK  11  j, 

6750  FOR  1=1  TO  LEN<A*)/2 
6760  PRINT  @D : "H" ; 

6770  NEXT  I 
6780  PRINT  @D:A$; 

6790  RETURN 

6800  PRINT  @D , 1 7 : 1 ,2,2 
6810  CHARSIZE  2 

6820  H$="Natural  Frequency  Hz" 

6830  MOVE  @D : < F7+F6 ) /2 , 0 
6840  PRINT  @D ; "JJJ"  ; 

6850  FOR  1=1  TO  LEN(H$)/2 
6860  PRINT  @D: "Hn \ 

6870  NEXT  I 
6880  PRINT  @D:H*> 

6890  GO  TO  R 1  OF  6900,6940,6960,6920,6980 
6900  Vi="Pseudovelocity  n/sec" 

6910  GO  TO  7010 

6920  V$  = "Rel .  Velocity  M/sec” 

6930  GO  TO  7010 

6940  V$="Abs.  Acceleration  g " 

6950  GO  TO  7010 

6960  V$=”Rel.  Displacement  mm11 
6970  GO  TO  7010 
6980  V*="  Phase  deg" 

6990  MOVE  @DiF6,9Q 

7000  GO  TO  7020 

7010  MOVE  @D : F6 , 0 . 5*1 NT  <  Ml ) 

7020  PRINT  @D: "HHHHHH” j 
7030  FOR  1=1  TO  LEN<V*)/2 
7040  PRINT  @D; ”KM ; 

7050  NEXT  I 

7060  FOR  1=1  TO  LEN(V$) 

7070  X$=SEG ( VS , I  ,  1  ) 

7080  PRINT  @D ; X$ ; ”HJ  "  j 
7090  NEXT  I 

7100  FOR  I=F6  TO  F7  STEP  <F7-F6>/10 

71 10  MOVE  @D: I >0 

7120  PRINT  <?D  i  "HJ  *  ;  INT  <  I  >  ; 

7130  NEXT  I 

7140  IF  R 1 =5  THEN  7270 
7150  IF  ABS ( Ml ) <  =  1 0  THEN  7190 
7160  IF  ABS (Ml ><-100  THEN  7210 
7170  IF  ABS(M1 ) <=1 000  THEN  7230 
7180  IF  ABS < Ml >>1000  THEN  7250 
7190  FOR  1=0  TO  ABS(Ml)  STEP  1 
7200  GO  TO  7280 

7210  FOR  1=0  TO  ABS(Ml)  STEP  10 
7220  GO  TO  7280 

7230  FOR  1=0  TO  ABS(Ml)  STEP  100 
7240  GO  TO  7280 

7250  FOR  1=0  TO  ABS<M1)  STEP  1000 


GO  TO  7280 

FOR  1  =  0  TO  180  STEP  45 
MOVE  @D : F6 , I 
PRINT  @D : "HHH " > I ; 

NEXT  I 

HOVE  @D:F6,0 
PRINT  8D , 1 7 : 1 ,2 
A*="K(Zero  Damping ) " 

PRINT  8D: M  J  J  J  J  "  ;  A$ 

1  =  1 

FOR  1*1  TO  CO 

IF  F5< I ) >0 . 01  THEN  7390 

NEXT  I 

MOVE  8D:F5(I) ,ABS(K(I) ) 

FOR  1  =  1  TO  CO 

IF  F5 ( I >  =  0 . Q 1  THEN  7360 

DRAW  @D : F5< I ) , ABS<K ( I ) ) 

NEXT  I 
RETURN 

REM  *****************  WRITE  DATA  *********************************> 
CHARSIZE  3 

PRINT  "  JLGENTER  FILE  NUMBER  (DEVICE  NR. 3)  TO  WRITE  DATA  ON:  - ; 

INPUT  D1 
FIND  83:D1 

WRITE  @3 : C 0 , W* , F5 , X5 , V5  , G5 , T5 , Y5 , Z5 
PRINT  ‘'RESULTS  STORED  IN  FILEG  *"  ;D1 
RETURN 

REM  ******************  READ  DATA  ********************************* 
CHARSIZE  3 

PRINT  "JLGENTER  FILE  NUMBER  (DEVICE  NR. 3)  TO  READ  DATA  FROM:  “ * 
INPUT  D3 

DELETE  F5,X5,V5,G5,T5,Y5,Z5,K,A 
FIND  @3:D3 
READ  83: CO 
FIND  83:D3 

DIM  F5(C0) , X5 ( CO ) , V5 ( CO ) , G5 ( CO ) , T5 ( CO ) , Y5 ( CO ) , Z5 ( CO ) ,  K (CO ) 

READ  83 : CO , W* , F5 , X5 , V5 > G5 , T5 , Y5 , Z5 

IF  W1 >0  THEN  7650 

PRINT  "DATA  ARE  IN  MEMORY ! G ■ 

RETURN 

REM  ***********************  LOG.  SPECTRA  ********************** 
VIEWPORT  5*1 .79,110-4*1 .79,10+3*2.82,100-3*2.82 
GO  TO  7740 

VIEWPORT  8*1 .79,130-4*1 .79,3*2.82,100-3*2.82 
IF  D=4  THEN  7720 
GO  TO  7740 

VIEWPORT  8*1 .79,150-4*1 .79,3*2.82,100-3*2.82 

PRINT  @D, 17: 1 ,2,1 .8 

CHARSIZE  1 

DELETE  K 

DIM  K ( CO ) 

CALL  "TIME" ,T$ 

PRINT  8D : "L">R* 

PRINT  @D:Ti; 

REM 

WINDOW  LGT ( N6) ,LGT(N7) , LGT ( N1 ) , LGT ( N2 ) 

IF  P  > 1 0  THEN  7980 
Hi =INT ( LGT  <  N1 ) ) 

FOR  1  =  1 0 AH1  TO  10A(H1  +  1 )  STEP  10AH1*P 
MOVE  8D : LGT ( N6 ) , LGT ( I > 


7860  DRAW  8D : LGT ( N7  > ,LGT < I ) 

7870  NEXT  I 
7880  HI =H1 +1 

7890  IF  HI < =LGT < N2 )  THEN  7840 
7900  HI =LGT <N6) 

7910  FOR  1  =  1 0  AH1  TO  10A<H1+1)  STEP  10AH1*P 
7920  MOVE  @D : LGT <  I )  , LGT ( N1 ) 

7930  DRAW  @D : LGT < I ) , LGT ( N2 ) 

7940  NEXT  I 
7950  HI  =H1 +1 

7960  IF  HI <  =LGT ( N7 )  THEN  7910 

7970  GO  TO  8020 

7980  AXIS  @D: t ,1 ,0,-2 

7990  MOVE  @D:LGT<N6) ,LGT(N2) 

8000  DRAW  @D : LGT (N7 ) ,  LGT ( N2) 

8010  DRAW  @D  s  LGT <N7) ,  LGT  <N1 ) 

8020  PRINT  @0,17:2,3 
8030  CHARSIZE  4 

8040  GO  TO  R  OF  8050,8070,8090 
8050  A$= "MAXIMAX  SHOCK  SPECTRUM" 

8060  GO  TO  8100 

8070  A*= "RESIDUAL  SHOCK  SPECTRUM" 

8080  GO  TO  8100 

8090  A$= "MAXIMAX  and  RESIDUAL  SHOCK  SPECTRUM" 

8100  GOSUB  8170 

8110  IF  W1 >0  THEN  8240 

8120  A$=W$ 

8130  PRINT  @D, 17: 1 .2,1 .8 
8140  CHARSIZE  2 
8150  GOSUB  8170 
8160  GO  TO  8240 

8170  MOVE  @D : LGT  <  1  0  A  <  0  . 5*LGT  <  N6*N7 ) ) ) , LGT ( N2  > 
8180  PRINT  @D : "KK"; 

8190  FOR  1=1  TO  LEN<A$>/2 
8200  PRINT  @D : "H"; 

8210  NEXT  I 
8220  PRINT  @D : AS ; 

8230  RETURN 

8240  GO  TO  R  OF  8250,8270,8290 
8250  V$="Pseudo velocity  n/sec" 

8260  GO  TO  8300 

8270  V$="Rel .  Velocity  «/secM 

8280  GO  TO  8300 

8290  V$= "Vel oc i t y  n/ sec" 

8300  MOVE  @D:LGT(N6) ,LGT( 10A < 0 .5*LGT(N2*N1 ) ) ) 

8310  PRINT  @0,17:1 .5,2.5 

8320  CHARSIZE  2 

8330  PRINT  @D: "HHHHHHH"; 

8340  FOR  1*1  TO  LEN<V*>/2 
8350  PRINT  @D : "K"; 

8360  NEXT  I 

8370  FOR  1=1  TO  LEN(V$) 

8380  X*=SEG<V*,I,1 ) 

8390  PRINT  @D:X$; "HJ"; 

8400  NEXT  I 

8410  MOVE  @D : LGT ( 1 0  A ( 0 . 5*LGT ( N6*N7 ) ) ) , LGT ( N1 ) 
8420  H$*"Natural  Frequency  Hz" 

8430  FOR  1*1  TO  LEN<H$>/2 
8440  PRINT  @D:  "H"  •, 

8450  NEXT  I 
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8460  PRINT  @D  :  "  J  J  J 11 ,  H$  j 
8470  HI -0 

8480  FOR  I=LGT<N1)  TO  LGT <N2 )  STEP  10AH1 
8490  MOVE  @D : LGT ( N6 )  ,  I 

8500  PRINT  ®D:  "HHHHHH;M10A,,;INT(H1+LGT<N1>) 

8510  HI  =H1  -**1 
8520  NEXT  I 
8530  HI =0 

8540  FOR  I =LGT ( N6 )  TO  LGT ( N7 )  STEP  10AH1 
8550  MOVE  3D: I,LGT<N1 ) 

8560  PRINT  @D: "HHJ" j "  1  0  A  M  ;  INT  <  HH-LGT  < N6 > ) 

8570  HI =H1 + 1 
8580  NEXT  I 

8590  MOVE  @D : LGT <  N6 )  , LGT < N1 ) 

8600  IF  D=32  THEN  8620 
8610  PRINT  @0,17:1,2 
8620  CHARSIZE  2 
8630  A$=” (Zero  Dawping)" 

8640  PRINT  @D : "JJJ" ;A$ 

8650  GOSUB  8770 
8660  GOSUB  8820 
8670  IF  R=3  THEN  8690 
8680  RETURN 
8690  R=R+1 

8700  IF  D<  >4  THEN  8740 

8710  PRINT  "LYOU  MAY  CHANGE  THE  PEN  COLOR  G • ; 

8720  PRINT  "  Press  RETURN  to  continue" 

8730  INPUT  AS 
8740  GOSUB  8770 
8750  GOSUB  8820 
8760  RETURN 

8770  GO  TO  R  OF  8780,8800,8780,8800 
8780  K  =V5 
8790  GO  TO  8810 
8800  K  = Y5 
8810  RETURN 
8820  1=1 

8830  FOR  1=1  TO  CO 

8840  IF  F5(I)>0.01  THEN  8870 

8850  NEXT  I 

8860  RETURN 

8870  MOVE  @D:LGT(F5(I)  ), LGT<ABS<K<D  )  ) 

8880  FOR  1=1  TO  CO 
8890  IF  F5 < I ) =0 . 01  THEN  8830 
8900  IF  R < *3  THEN  8920 
8910  DASH  170 

8920  DRAW  9D : LGT <F5< I ) ) ,LGT ( ABS(K < I )  )  ) 

8930  NEXT  I 

8940  IF  W1 >0  THEN  8960 
8950  DASH  0 
8960  HOME 
8970  RETURN 

8980  REM  ******»#*»*»**  NUMERICAL  OUTPUT  ##***#***#****#*«#*«**#»****** 

8990  WINDOW  0,130,0,100 

9000  VIEWPORT  0,130,0,100 

9010  11*0 

9020  12*1 

9030  N* 1 

9040  IF  D*32  THEN  9070 
9050  IF  D*40  THEN  9130 
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9060  PRINT  90,1711 .5,2.5 
9070  CHARSIZE  3 
9080  PRINT  80:  “l," 

9090  HOVE  80:65,99 
9100  FOR  1=1  TO  LEN<US>/2 
9110  PRINT  80: "H"; 

9120  NEXT  I 
9130  PRINT  80 : U* 

9140  CHARSIZE  1 
9150  CALL  "TIME" ,T$ 

9160  PRINT  80:  "IT"  jR$, 'V>T*,  **JIM 
9170  IF  0=40  THEN  9190 

9180  PRINT  80:  USING  9190  :  "KKK  H  ;  "PAGE**  ;  12  j  11 J  J" 

9190  IHAGE  3 A , 1 22X , 4A , 20 , 2A 
9200  CHARSIZE  2 

9210  PRINT  80 :  ”  FREQ  SHOCK  SPECTRUM  FOURIER  TRANS** 

9220  PRINT  80:  “  <HZ)  < HM)  <H/S>  (G)  <MS)  <M/S)  (DEG)^** 

9230  IMAGE  40 . 1 0 , X , 50 . 20 , X , 40 , 20 , X , 5D . 1 0 , X , 40 . 20 , X , 4D . 20 , 2X , 40 
9240  IMAGE  65T , 40 . 1 D , X , 50 . 20 , X , 40 . 20 , X , 50 . 1 D , X , 40 . 20 , X , 40 . 20 , 2X , 4D 
9250  FOR  1=1+11  TO  CO 
9260  IF  I >50*N  THEN  9300 

9270  PRINT  80:  USING  9230 : F5< I ) ,X5< I ) , V5< I > ,G5( I ) ,T5< I ) , Y5( I ) , Z5< I ) 

9280  NEXT  I 
9290  RETURN 

9300  IF  0=40  THEN  9350 

9310  CHARSIZE  1 

9320  PRINT  80: MA JJJJJ" 

9330  CHARSIZE  2 
9340  PRINT  80: "J" 

9350  FOR  1*1 +50*N  TO  CO 
9360  IF  I > ( N+l ) *50  THEN  9430 
9370  IF  0=40  THEN  9400 

9380  PRINT  80:  USING  9240 :F5< I ) , X5< I ) , V5( I > ,G5( I) ,T5( I > , Y5( I ) ,Z5( I ) 

9390  GO  TO  9410 

9400  PRINT  80:  USING  9230 : F5  < I > , X5 < I > , V5 ( I > , G5 ( I > , T5 ( I ) , Y5< I > , Z5 ( I ) 

9410  NEXT  I 

9420  RETURN 

9430  IF  0*40  THEN  9450 

9440  CALL  HUAITH,4 

9450  11*100 

9460  12*12+1 

9470  N-N+2 

9480  IF  0*40  THEN  9500 
9490  GO  TO  9040 
9500  GO  TO  9250 

9510  REM  *••»»*••***«•«*•****»  ANALYSIS  ******************************** 

9520  R-3 

9530  GOSUB  7670 

9540  R*1 

9550  MOVE  80:3.11,2 

9560  PRINT  80;  USING  "2  A , 5X , 4A , 3X , 7A **  :  "NR  "  ,  C HZ 3  "  ,  **  L M/SEC  1  H 
9570  POINTER  X,Y,Q* 

9580  IF  Q$* "0  H  THEN  9750 

9590  MOVE  80 : X-LGT ( 1 , 022 > , Y+LGT ( 1 , 022 ) 

9600  CHARSIZE  2 

9610  PRINT  80:  ,,JA**;R 

9620  GOSUB  9650 

9630  R-R+l 

9640  GO  TO  9570 

9650  VIEWPORT  106,130,0,100 


42 


9660  UINDOU  0,32,0,100 
9670  MOVE  @0:0,94 
9680  FOR  1*1  TO  R 
9690  PRINT 
9700  NEXT  I 

9710  PRINT  80s  USING  9720 : " J A " ,R , " : " , 1 0 AX ,  1  0 A Y 
9720  IMAGE  2A , 2D , A , 2X , 4D . 1 D , 3X , 3D . D 

9730  VIEWPORT  5*1.79,110-4*1.79,10+3*2.82,100-3*2.82 
9740  WINDOW  LGT ( N6 ) ,LGT < N7 ) ,LGT <N1 ) , LGT <N2) 

9750  DASH  0 
9760  HOME 
9770  RETURN 

9780  REM  *********************  MULTI  DISPLAY  *************************** 
9790  CHARSIZE  3 
9800  D-32 

9810  PRINT  "LfENTER  NUMBER  OF  DESIRED  SPECTRAG; 

9820  INPUT  Wl 
9830  DIM  W2 ( Wl ) 

9840  PRINT  ".ENTER  FILE  NUMBERS  ON  DEVICE  NR  .  3G  : 

9850  INPUT  U2 

9860  PRINT  ".ENTER:.  1  FOR  MAXIMAX-SPECTRA. 2  FOR  RES I DUAL -SPECTRAG :  “  > 
9870  INPUT  R 
9880  D3-W2  < 1 ) 

9890  COSUB  7570 

9900  GOSUB  7690 

9910  GOSUB  10060 

9920  IF  Wl-1  THEN  10050 

9930  FOR  10-2  TO  Wl 

9940  D3=W2( 10 ) 

9950  GOSUB  7570 
9960  DASH  11*10 
9970  GOSUB  8770 
9980  GOSUB  8820 
9990  GOSUB  10060 
10000  NEXT  10 
10010  DASH  0 
10020  Wl -0 

10030  DELETE  W2 , X , Y , 1 0 , D3 , Q* 

10040  W*- "TEST " 

10050  RETURN 

10060  POINTER  X,Y,Q* 

10070  MOVE  9D : X , Y 
10080  DRAW  @D : X+0 . 25 , Y 
10090  PRINT  9D : "  " jW$ 

10100  RETURN 
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1  REM  /////////////////  8SPRGM/PART  III  \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\ 

2  INIT 

3  GO  TO  260 

4  IF  0*3  THEN  94 

5  GO  TO  7450 

8  IF  0=3  THEN  94 

9  GO  TO  7530 

12  IF  0=3  THEN  94 

13  GO  TO  8980 
16  GO  TO  190 
20  GO  TO  100 

24  IF  0=3  THEN  94 

25  R  =  1 

26  GO  TO  7690 

28  IF  0=3  THEN  94 

29  R«2 

30  GO  TO  7690 

32  IF  0=3  THEN  94 

33  R*3 

34  GO  TO  7690 

36  IF  0=3  THEN  94 

37  GO  TO  9510 

40  IF  0=3  THEN  94 

41  CO  TO  9780 
45  RETURN 

48  IF  0=2  THEN  90 

49  GO  TO  510 

50  DELETE  A 

52  IF  0=2  THEN  90 

53  GO  TO  1780 

54  DELETE  A 

64  IF  0=3  THEN  94 

65  R 1  =  1 

66  GO  TO  6010 

68  IF  0=3  THEN  94 

69  R 1 =4 

70  GO  TO  6010 

72  IF  0=3  THEN  94 

73  R 1 =2 

74  GO  TO  6010 

76  IF  0=3  THEN  94 

77  R 1 =3 

78  GO  TO  6010 

80  IF  0=3  THEN  94 

81  R 1 =5 

82  GO  TO  6010 

90  PRINT  "You  got  to  PRESS  KEY  5  first»G" 

92  RETURN 

94  PRINT  "You  got  to  PRESS  KEY  4  firstIG" 

96  RETURN 

100  REM****************  LOAD  ROUTINE  *********************************** 
110  IF  0=3  THEN  170 

120  DELETE  F5 , X5 , V5 , G5 , T5 , Y5 , Z5 , C ,C0 , K  ,  A 
130  C0=90 
140  0=3 

150  APPEND  "8SPRGM/P ARTI I I  *  ;  1 

160  PRINT  " JG5PRGM/PART  III  IS  LOADED ! " 

170  DELETE  5140,20000 
180  RETURN 
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190  IF  0-2  THEN  240 
200  0*2 

210  DELETE  270,20000 

220  APPEND  "8SPRGM/P ARTI I "  ;260 

230  RENUMBER  6000,10,260 

240  PRINT  M JGSPRGM/P ART  II  IS  LOADED!'1 

250  RETURN 

260  REM**************  DEFAULT  CONDITION  ******************************** 
270  WS= "TEST " 

280  0-3 
290  Wl-0 
300  D*32 
310  Fl*l 
320  F2*l 000 
330  F8*l 000 
340  N1 *0 . 0 1 
350  N2*l 0 0 
360  N6*l 
370  N7«l 000 
380  C0*90 
390  N0*4096 
400  K 1*2 
411  K  2*2 
420  K3*l 0 
430  P=1 0 

440  R$*"NRL  -  CODE  5842" 

450  DIM  M9(7,2) , A2  <  7 , 2 ) , C ( C 0 ) 

460  M9=0 
470  A2*0 

480  DIM  Df  (8) , Ht<20 ) , V$<22) ,G$<7) ,M$(7> , Nf<8> ,X$(22) , P*(35) ,L$(4) ,  Q$<  1 ) 
490  PRINT  ”  GSPRGM/PART  III  IS  LOADED  (DEFAULT  CONDITIONS)!" 

500  RETURN 

510  REM  *****************  READ  RECORD*********************************** 
520  CHARSIZE  3 

530  PRINT  "LJENTER  FILE  NUMBER  (ON  DEVICE  3)  TO  READ  RECORD  FROMG ; 

540  INPUT  D1 
550  DIM  A < NO ) 

560  W3*130 

570  IF  D<  >4  THEN  600 
580  W3*150 

590  PRINT  8D , 1 7 : 1 .2,1.8 
600  CHARSIZE  2 
610  CALL  "TIME "  ,  T$ 

620  PRINT  8D:  m\xK"  ,R$ 

630  PRINT  8D :  T$ 

640  FIND  83 : D1 

650  READ  83 : W$ , R2 , H , T 1 , M 

660  FIND  83 ; D1 

670  READ  83 : W$ , R2 , H , T1 ,M,M9(M,2) ,A2(M,2) ,A 
680  M9<M,1)-M 
690  A2(M,1)»M 
700  S*M9 < M , 2 > 

710  WINDOW  0,W3, 0,100 
720  VIEWPORT  0,W3, 0,100 
730  IF  W1>0  THEN  810 
740  PRINT  8D, 17:2,3 
750  CHARSIZE  4 
760  MOVE  8D ; W3/2 , 98 
770  FOR  I»1  TO  LEN<Wi)/2 
780  PRINT  8D : "H"  j 
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790  NEXT  I 
800  PRINT  @D:W* 

810  PRINT  @0,17:1 .2,1.8 
820  CHARSIZE  2 

830  GO  TO  R2  OF  840,970,1110 

840  GQSUB  1280 

850  S1=S 

860  GOSUB  1460 

870  PRINT  @0 : "# " 

880  CALL  "DIFS", A, A 

890  S1=S1/H 

900  GOSUB  1340 

910  GOSUB  1460 

920  CALL  "DIF3  “  ,  A  ,  A 

930  SI =S1  *1 QA3/(H*9,8044) 

940  GOSUB  1400 
950  GOSUB  1460 
960  GO  TO  1230 
970  GOSUB  1340 
980  S1=S 
990  GOSUB  1460 
1  000  PRINT  @D  : 

1010  CALL  MINT" , A, A 

1020  S 1 *S1 *H 

1030  GOSUB  1280 

1040  GOSUB  1460 

1050  CALL  "DIF3 " , A , A 

1060  CALL  M  D I F  3  "  ,A,A 

1070  S1=*S1*10A3/(HA2*9.8044> 

1080  GOSUB  1400 
1090  GQSUB  1460 
1100  GO  TO  1230 
1110  GOSUB  1400 
1120  S1=*S 
1130  GOSUB  1460 
1140  PRINT  @D: "a" 

1150  CALL  "INT "  ,  A  ,  A 

1160  S1=S1*1 . QE-3*H*9.8044 

1170  GOSUB  1340 

1180  GOSUB  1460 

1190  CALL  "INT  "  ,  A  ,  A 

1200  SI =S1 #H 

1210  GOSUB  1280 

1220  GOSUB  1460 

1230  VIEWPORT  0,W3, 0,100 

1240  WINDOW  0,W3, 0,100 

1250  MOVE  CD: 0,0 

1260  PRINT  @D:M*  orig.  ti«e  record”; 
1270  RETURN 

1200  VIEWPORT  5 , W3-30 ,4,32 
1290  V$» "DISPLACEMENT " 

1300  G**"D  wax  : " 

1310  M$-"D  min  : " 

1320  Nf * w  [**]" 

1330  RETURN 

1340  VIEWPORT  5 , W3-30 , 36 , 64 
1350  V$* "VELOCITY " 

1360  G**"V  nax 
1370  M$*"V  «in  :  " 

1380  im/%mc J H 
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1390  RETURN 

1400  VIEWPORT  5,W3-30,68,96 
1410  VS» "ACCELERATION 11 
1420  GS»”A  Max  : 

1430  MS* "A  Min  : “ 

1440  NS* "  C q ] " 

1450  RETURN 

1460  CALL  "MAX " , A , VI ,14 
1470  CALL  "MIN” , A , V2 > 13 
1480  WINDOW  0 , NO , V2 > VI 
1490  IF  D*32  THEN  1550 
1500  MOVE  @D : 0 , 0 
1510  FOR  1*1  TO  NO 
1520  DRAW  @D : I , A ( I ) 

1530  NEXT  I 

1540  GO  TO  1560 

1550  CALL  "D ISP  ",  A 

1560  AXIS  @D : 40  0 , 0 

1570  FOR  1*0  TO  NO  STEP  2000 

1580  MOVE  @D : 1 , 0 

1590  PRINT  @D : "HJ",I*H; 

1600  NEXT  I 

1610  PRINT  @D ; "HHHHHHHHHHH CMS]" 

1620  MOVE  @D : 420 0 , VI 

1630  PRINT  @D:  USING  1670: " J " ,GS  > VI #S1 , NS ; 

1640  MOVE  @D: 4200, VI 
1650  PRINT  @D: "J” j 

1660  PRINT  8D:  USING  1 670 : " J " , MS , V2*S1 , NS ; 

1670  IMAGE  A , 7A , 6D . 1 D ,  1 1 A 
1680  MOVE  @D ; 0 , 0 . 5*  < V2+V1 ) 

1690  PRINT  6D: ”HH " ; 

1700  FOR  1*1  TO  LEN<VS)/2 
1710  PRINT  @D : "K 
1720  NEXT  I 

1730  FOR  1*1  TO  LEN(VS) 

1740  XS»SEG <  VS , 1 , 1  ) 

1750  PRINT  @D:XS; ”HJ"> 

1760  NEXT  I 
1770  RETURN 

1780  REM  ********************  PROCESS  ********* ******* * ****** *********** 
1790  CHARSIZE  3 
1800  CALL  "TIME ”, BS 
1810  PRINT  9D: ”t";BS; "i" 

1820  PRINT  @32,26:2 
1830  IF  FI < 1  THEN  1890 
1840  IF  FI >1000  THEN  1890 
1850  IF  F2< 1  THEN  1900 
1860  IF  F2 > 1 000  THEN  1900 
1870  IF  F2<  *F 1  THEN  1900 
1880  GO  TO  1930 
1890  IF  C0=1  THEN  1930 

1900  PRI  "ANALYSIS  FREQUENCY  IS  OUT  OF  RANGE  ENTER  NEW  VALUES  OF  F1,F2G" 
1910  INPUT  F1,F2 
1920  GO  TO  1800 

1930  DELETE  F5 , X5 , V5 , G5 , T5 , Y5 , Z5 , C , K , A 
1940  DIM  A< NO ) 

1950  C3*INT < (MEMORY/8-1 8 ) /8 ) 

1960  DIM  F5<C0> ,X5<C0> ,V5(C0) ,G5<C0) ,T5<C0) ,Y5(C0> ,Z5(C0> ,C(C0) 

1970  FIND  63 : D1 

1980  READ  83:WS,R2,H,T1 ,M,S 
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1 990  FIND  @3 ! D1 

2000  E1=INT(1000/(K1*F8*H>> 

2010  E2=INT< 1000/(K2*F8*H) > 

2020  E3=INT(1000/(K3*F8*H> > 

2030  F5( 1 ) =F 1 

2040  IF  C0  =  1  THEN  2080 

2050  FOR  1=2  TO  CO 

2060  F5( I )=(F2/F1 )*<1/<C0-1 > )*F5< 1-1 > 

2070  NEXT  I 
2080  X5=0 . 0 1 
2090  V5=0 . 0 1 
2100  G5=0 . 0 1 
2110  T5=0 , 0 1 
2120  Y5=0 . 0 1 
2130  Z5=0 . 0 1 
2140  11=1 
2150  12=0 

2160  IF  FI < 1 0  THEN  2200 
2170  12=1 

2180  IF  FI < 1 00  THEN  2200 

2190  12=2 

2200  GO  TO  2370 

2210  E=E1 

2220  GOSUB  2460 

2230  IF  I1>C0  THEN  2310 

2240  GO  TO  2370 

Torn  tr-rrp 

2260  GOSUB  2460 

2270  IF  1 1  )C0  THEN  2310 

2280  GO  TO  2370 

2290  E=E3 

2300  GOSUB  2460 

2310  DELETE  A 

2320  CALL  "TIME  ", T$ 

2330  GOSUB  3320 
2340  PRINT  @32,26:0 

2350  PRINT  @D:"_MAX.  POSSIBLE  ANALYSIS  FREQUENCIES:  ;C3 

2360  RETURN 
2370  F3=F5< I 1 ) 

2380  C1=I1 

2390  12=12+1 

2400  FOR  11=1  TO  CO 

2410  IF  F5(I1)>10AI2  THEN  2430 

2420  NEXT  II 

2430  F4=F5< 1 1  —  1  > 

2440  GO  TO  12  OF  2210,2250,2290 
2450  RETURN 

2460  READ  @3 : U* , R2 , H , T1 , M , S , A2<M , 2 > , A 

2470  FIND  @3:D1 

2480  IF  E<  =1  THEN  2550 

2490  J=1 

2500  FOR  1=1  TO  NO  STEP  E 
2510  A<  J  >=A< I > 

2520  J=J+1 
2530  NEXT  I 
2540  GO  TO  2560 
2550  E=1 
2560  N5*N0/E 
2570  H2=H»E 

2580  GO  TO  R2  OF  2590,2630,2660 
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2590  CALL  ,,DIF3,,>A>A 

2600  CALL  l'DIF3"  ,  A ,  A 

2610  S2=S 

2620  GO  TO  2670 

2630  CALL  WDIF3",A>A 

2640  S2=S*H2 

2650  GO  TO  2670 

2660  S2=S*10A-3*H2A2*9.8044 

2670  C2=I 1-1 

2680  Pl  =  0 . 0062831853*F3*H2 
2690  FOR  I=C1  TO  C2 
2700  C ( I ) =2*C0S (PI) 

2710  IF  Cl =C2  THEN  2730 

2720  P1=P1*(F4/F3)A(1/<C2“C1 )  ) 

2730  NEXT  I 

2740  IF  P 1 <  0 . 902  THEN  2790 
2750  IF  1 43 . 566/H2 >F2  THEN  2790 

2760  PRINT  "LGWARNING :  FREQUENCIES  ARE  TOO  HIGH  FOR  YOUR  TIME  INCREMENT" 
2770  PRINT  ”  ERRORS  ARE  OUTSIDE  THE  RANGE  -7  TO  +3  PERCENT  FOR H 

2780  PRINT  •  FREQUENCIES  ABOVE  " )  1 43 . 566/H2  ;  "  HERTZ" 

2790  FOR  I=C1  TO  C2 
2800  F=F5 ( I ) 

2810  IF  F <  =  1  THEN  2840 
2820  N=1000/(H2*5*LGT (F) ) 

2830  IF  N< N5  THEN  2850 

2840  N=N5 

2850  X1=0 

2860  X2=0 

2870  R 3=0 

2880  R 4=0 

2890  IF  F=0  THEN  2920 
2900  N4=N-1-50Q/<F*H2> 

2910  IF  N4  >3  THEN  2930 
2920  N4=3 

2930  FOR  J*1  TO  N-l 
2940  X3=C(I)*X2-XH-A(  J) 

2950  IF  J < N4  THEN  3000 

2960  IF  ABS ( X3 ) <  ABS ( R 4 )  THEN  3030 

2970  R4=*X3 

2980  T 4= J 

2990  GO  TO  3030 

3000  IF  ABS(X3XABS<R3)  THEN  3030 

3010  R  3=X3 

3020  T3= J 

3030  X 1 =X2 

3040  X2*X3 

3050  NEXT  J 

3060  IF  ABS < R3 ) > ABS ( R4 )  THEN  3090 

3070  R3=R4 

3080  T3=T 4 

3090  X=S2*R3 

3100  V=0 . 00628318*F*X 

3110  G=6 . 28831 8*F*V/9 ,8044 

3120  T=T1 +(T3-0 , 5>*H2 

3130  Y  =  0  . 006283 1 8*F*S2*R4 

3140  Z=1  . 0E-3* ( T 1  +  <  T4-0 , 5 ) *H2 > *F  +  0 , 25 

3150  Z=Z- INT < Z ) 

3160  IF  Y  >  0  THEN  3180 
3170  Z=Z-0 . 5 

3180  IF  Z<*0 .5  THEN  3200 
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3190  Z=Z-1 
3200  Z=360*Z 
3210  F5<I)=F 
3220  X5<I)=X 
3230  U5(I)=V 
3240  G5  < I )*G 
3250  T5 ( I >  =T 
3260  Y5(I)=Y 
3270  Z5(I)=Z 

3280  PRINT  @D:  USING  3290 : I ,F,U,N; 

3290  IMAGE  20T, 3D, 4X , 4D . ID, 4X, 5D . 2D.3X.6D 
3300  NEXT  I 
3310  RETURN 

3320  REM  ******  SUBROUTINE  TIME  ************************************** 
3330  DIM  B( 2 ) 

3340  C$=B* 

3350  FOR  1=1  TO  2 
3360  X*=5EG(C$,11 ,8) 

3370  Y*=SEG(X$,4,2> 

3380  Z*=SEG<X$,7,2) 

3390  B< I >=UAL(X$>+( VAL( Y$>*60+VAL(Z*>  >/3600 
3400  C$*T $ 

3410  NEXT  I 
3420  U=B(2)-B(1) 

3430  B3=» < U-INT < U )  )*60 
3440  B4=< B3-INT < B3 ) ) *60 

3450  PRINT  8D : “ _ CALCULATION  TIME  ( hr* : win : sec ) : 

3460  PRINT  @D:  USING  3470  s INT <U INT < 83 ), M , INT ( B4 ) 

3470  IMAGE  2D, A, 2D, A, 2D 
3480  RETURN 

5000  REM  - 

5010  REM - - 

5020  REM  THIS  PROCEDURE  MAY  BE  USED  TO  CALCULATE  SHOCK  SPECTRA  OF 

5030  REM  TIME  RECORDS  STORED  ON  EXT.  TAPE  (FILE  1  TO  7)  AND  TO  STORE 

5040  REM  THE  ACCORDING  DATA  ON  INT.  TAPE  (FILE  *  TO  FILE  *+7>. 

5050  FOR  DIM  TO  7 
5060  GOSUB  52 
5070  N=Dl+0 
5080  FIND  N 

5090  MARK  1 ,7000 

5100  FIND  N 

5110  WRITE  C0,W*,F5,X5,V5,G5,T5,Y5,Z5 
5120  NEXT  D1 
5130  RETURN 
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